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Chapter 1 
Introduction 



Problems of interface growth have received much attention recently [1-3]. Such are, for 
example, the duffusion limited aggregation (DLA) [4], random sequential adsorption (RSA) 
[5], Laplacian growth [6-8] or flame front propagation [9]. We will mainly pay attention in 
this Thesis to the numerical and analytical investigation of the last two problems. In addition 
to the fact that flame front propagation is an interesting physical problem we feel that we 
can also explain experimental results on the basis of theoretical investigations. There exists 
possibility to use methods found for the flame front propagation, in different fields where 
similar problems appear such as the important model of Laplacian growth . 

The premixed flame - the self-sustaining wave of an exothermic chemical reaction - is 
one of the basic manifestations of gaseous combustion. It is well established, however, that 
the simplest imaginable flame configuration - unbounded planar flame freely propagating 
through initially motionless homogeneous combustible mixture - is intrinsically unstable and 
spontaneously assumes a characteristic two- or three-dimensional structure. 

In the recent paper of Gostintsev, Istratov and Shulenin [10] an interesting survey of 
experimental studies on outward propagating spherical and cylindrical flames in the regime 
of well developed hydrodynamic (Darrieus-Landau) instability is presented. The available 
data clearly indicate that freely expanding wrinkled flames possess two intrinsic features: 

1. Multi- quasi- cusps structure of the flame front. (The flame front consists of a large 
number of quasi-cusps, i.e., cusps with rounded tips.) 

2. Noticeable acceleration of the flame front 

Moreover, the temporal dependence of the flame radius is nearly identical for all premix- 
tures discussed and correlates well with the simple relation: 



Here -Ro(^) is the effective (average) radius of the wrinkled ffame and A, B are empirical 
constants. 

In this Thesis we study the spatial and temporal behavior of a nonlinear continuum model 
(i.e., a model which possesses an inflnite number of degrees of freedom) which embodies 
all the characteristics deemed essential to premixed flame systems; namely, dispersiveness, 
nonlinearity and linear instability. Sivashinsky, Filyand and Prankel [11] recently obtained 



Ro(t) = At^/^ + B 
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an equation, denoted by SFF in what follows, to describe how two-dimensional wrinkles of 
the cylindrical premixed flame grow as a consequence of the well-known Landau-Darrieus 
hydrodynamic instability. The SFF equation reads as follows: 



dR _ U, dRV _DM_d^ ^rim + u n 2) 

dt ~ 2R^\t) [dej ^ Ro\t) de^ + 2it:oW^^ ' ^ ^ ^ 

where 0< ^ < 27r is an angle, R(^,t) is the modulus of the radius- vector on the flame interface, 
Ub, Dm, 7 are constants. 

1 /•27r 

UR) = -y^n cos[n{9 - 9*)]R{9*, t)d9* = 



R<,(t) = ^ r R(»,t)de . (1.4) 

Sivashinsky, Filyand and Frankel [11] made a direct numerical simulation of this nonlinear 
evolution equation for the cylindrical flame interface dynamics. The result obtained shows 
that the two mentioned experimental effects take place. Moreover, the evaluated acceleration 
rate is not incompatible with the power law given by eq.(l.l). For comparison, numerical 
simulations of freely expanding diffusively unstable flames were presented as well. In this 
case no tendency towards acceleration has been observed. 

In the absence of surface tension, whose effect is to stabilize the short-wavelength pertur- 
bations of the interface, the problem of 2D Laplacian growth is described as follows 

(a^ + (?> = o. (1.5) 

u |r(t)= ,dnU |s= 1 • (1.6) 

Vn = dnU |r(t) . (1.7) 

Here u{x,y;t) is the scalar field mentioned, T{t) is the moving interface, S is a fixed 
external boundary, dn is a component of the gradient normal to the boundary (i.e. the 
normal derivative), and Vn is a normal component of the velocity of the front. 

To obtain results for radial flame growth it is necessary to investigate the channel case 
first. The channel version of equation for flame front propagation is the so-called Michelson- 
Sivashinsky equation [12,13] and looks like 

I{E) = — P / dx* . (1.9) 

TT J —CO X — X 
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Figure 1.1: Giant cusp solution 

with periodic boundary condition on the interval x [0,L], where L is size of the system, u 
is constant, v > H is the height of the flame front point, P J is the usual principal value 
integral. 

Equations for flame front propagation and Laplacian growth with zero surface tension 
have remarkable property : these equations can be solved in terms of poles in the complex 
plane [6, 12, 14-16]. So we obtain a set of ordinary differential equations for the coordinates 
of these poles. The number of the poles is constant value in the system, but to explain such 
effect as growth of the velocity fiame front we need to consider some noise that is a source 
of new poles. So we need to solve the problem of interaction of the random fluctuations and 
the pole motion. 

The simplest case is the channel geometry. Main results for this case is existence of the 
giant cusp solution [12](Fig.l.l), which is represented in configuration space by poles which 
are organized on a line parallel to the imaginary axis. This pole solution is an attractor for 
pole dynamics. 

A complete analysis of this steady-state solution was first presented in Ref. [12] and the 
main results are summarized as follows: 

1. There is only one stable stationary solution which is geometrically represented by a giant 
cusp (or equivalently one finger) and analytically by N{L) poles which are aligned on 
one line parallel to the imaginary axis. The existence of this solution is made clearer 
with the following remarks. 

2. There exists an attraction between the poles along the real line. The resulting dynamics 
merges all the x positions of poles whose y-position remains finite. 

3. The y positions are distinct, and the poles are aligned above each other in positions 
Uj^i < Hj < i/j+i with the maximal being yN{L)- This can be understood from equations 
for the poles motion in which the interaction is seen to be repulsive at short ranges, 
but changes sign at longer ranges. 
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If one adds an additional pole to such a solution, this pole (or another) will be pushed to 
infinity along the imaginary axis. If the system has less than N{L) poles it is unstable 
to the addition of poles, and any noise will drive the system towards this unique state. 
The number N{L) is 

N{L)=\U- + l)] , (1.10) 



where 
with po 



is the integer part and 2ttL is a system size. To see this consider a system 

es and such that all the values of yj satisfy the condition < yj < y,nax- 
Add now one additional pole whose coordinates are Za = {xa, ya) with ya ^ ymax- From 
the equation of motion for ya, we see that the terms in the sum are all of the order of 
unity as is also the cot(ya) term. Thus the equation of motion of ya is approximately 

dya ^ 2N+1 1 

The fate of this pole depends on the number of other poles. If is too large the pole 
will run to infinity, whereas if is small the pole will be attracted towards the real 
axis. The condition for moving away to infinity is that N > N{L) where N{L) is given 
by (1.10). On the other hand the y coordinate of the poles cannot hit zero. Zero is a 
repulsive line, and poles are pushed away from zero with infinite velocity. To see this 
consider a pole whose yj approaches zero. For any finite L the term coth(|/j) grows 
unboundedly whereas all the other terms in the equation for the poles motion remain 
bounded. 

5. The height of the cusp is proportional to L. The distribution of positions of the poles 
along the fine of constant x was worked out in [12]. 

We will refer to the solution with all these properties as the Thual-Frisch-Henon (TFH)-cusp 
solution. 

The main results of our own work are as follow. Traditional linear analysis was made 
for this giant cusp solution. This analysis demonstrates the existence of negative eigenvalues 
that go to zero when the system size goes to infinity. 

1. There exists an obvious Goldstone or translational mode with eigenvalue Aq = 0. This 
eigenmode stems from the Galilean invariance of the equation of motion. 

2. The rescaled eigenvalues (L^Aj) oscillate periodically between values that are L-independent 
in this presentation. In other words, up to the oscillatory behavior the eigenvalues de- 
pend on L like L~^. 

3. The eigenvalues Ai and A2 hit zero periodically. The functional dependence in this 
presentation appears almost piecewise linear. 

4. The higher eigenvalues also exhibit similar qualitative behaviour, but without reaching 
zero. We note that the solution becomes marginally stable for every value of L for which 
the eigenvalues Ai and A2 hit zero. The dependence of the spectrum indicates that 
the solution becomes more and more sensitive to noise as L increases. 
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It was proved that arbitrary initial conditions can be written in the term of poles in the 
complex plane. Inverse cascade process of giant cusp formation was investigated numerically 
and analytically. Dependences of the flame front width and mean velocity were found. The 
next step in investigation of the channel case was the influence of random noise on the pole 
dynamics. The main effect of the external noise is the appearance of new poles in the minima 
of the flame front and the merging these poles with the giant cusp. The dependence of the 
mean flame front velocity on the noise and the system size was found. The velocity is almost 
independent on the noise until the noise achieves some critical value. In the dependence of 
the velocity on the system size we see growth of the velocity with some exponent until the 
velocity achieves some saturation value. 

Denoting v as the velocity of the flame front and L the system size: 

1. We can see two different regimes of behavior the average velocity v as a, function of 
noise / for flxed system size L. For the noise / smaller then same flxed val ue fcr 

Vr^f^. (1.12) 

For these values of / this dependence is very weak, and ^ 0.02. For large values of / 
the dependence is much stronger 

2. We can see growth of the average velocity v as a function of the system size L. After 
some values of L we can see saturation of the velocity. For regime / < fcr the growth 
of the velocity can be written as 

vr^L", // 0.35 ± 0.03 . (1.13) 

The dependence of the number of poles in the system and the number of the poles that 
appear in the system in unit time was investigated numerically as a function of the noise and 
the system parameters. The hfe time of a pole was found numerically. Theoretical discussion 
of the effect of noise on the pole dynamics and mean velocity was made [17]. 

Pole dynamics can be used also to analyse small perturbation of the flame front and make 
the full stability analysis of the giant cusp. Two kinds of modes were found. The flrst one 
is eigenoscillations of the poles in the giant cusp. The second one is modes connected to the 
appearance of the new poles in the system. The eigenvalues of these modes were found. The 
results are in good agreement wit h the traditional stability analysis [18]. 

The results found for the cannel case can be used to analyse flame front propagation in 
the radial case [19,20]. Main feature of this case is a competition between attraction of the 
poles and expanding of the flame front. So in this case we obtain not only one giant cusp but 
a set of cusps. New poles that appear in the system because of the noise form these cusps. 
On the basis of the equation of poles motion we can flnd connection between acceleration of 
the flame front and the width of the interface. On the basis of the result for mean velocity in 
the channel case the acceleration of the flame front can be found. So we obtain full picture 
of the flame front propagation in the radial case. 

The next step in the investigation of the problem is considering Laplacian growth with 
zero surface tension that also has pole solutions. In the case of Laplacian growth we obtain 
result that is analogous to the merging of the poles in the channel case of the flame front 
propagation: all poles coalesce into one pole in the case of periodic boundary condition or 
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two poles on the boundaries in the case of no-flux boundary conditions. This result can be 
proved theoretically [21]. 

In papaers [22-26] self-acceleration without involvement of the external forcing is con- 
sidered. No self- acceleration exist for the finite number of poles. So we can explain the 
self- acceleration and the appearance of new poles or by the noise or by the "rain" of poles 
from the "cloud" in infinity. Indeed, any given initial condition can be written as a sum of 
infinite number of poles( Sec. 2.4.1). Let us consider one pole that appears from the "cloud" 
in infinity. We neglect by the repulse force from the rest of poles in the system and consider 
only attraction force in Eqs.(4.13) —-^ For tq we can write in the case of self- acceleration 
'^o(''') — {a + tY, '^o(O) — , (3 > 1. So from r = to r = oo pole comes down a distance 

— Io° 2r]{T) ^'^ ~ — ^T^- So the "rain" come down the finite distance after the infi- 

r-o{0) f 

nite time and this distance converges to zero if ro(0) i— > oo! So we think that the appearance 

of new poles from the infinity can be explained only by the external noise. The characteristic 

_i 

size of cusp in the system £ ~ rn . So from Fig 2.23 the noise / ~ i ~ ^ is necessary 

for the appearance of new cusps in the system. If the noise is larger than this value the 
dependence on the noise is very slow( /°.2 for regime II and /°.02 for regime III). This result 
explains the weak dependence of numerical simulations on the noise reduction ( [22], Fig.2). 

Joulin et al. [27-30] use a vary similar approach for the channel and radial flame growth. 
But the main attention in our work is made to the velocity of flame front (self-acceleration for 
the radial case) and the flame front width. Main attention in the channel case in Joulin's work 
was made to the investigation of mean-spacing between cusps (crests). For the radial case 
only the linear dependence of radius on time (no self-acceleration) is considered in Joulin's 
work. Our works greatly complement each other but don't compete with each other. For 
example, for the distance between cusps (the mean value of cusp) without any proof we use 
Eq.(2.70). Fig.9( [28]) give us a excellent proof of this equation. 

The structure of this Thesis is as follow. Chapter 1 is this Introduction. 

In Chapter 2 we obtain main results for the channel case of the flame front propagation. 
We give results about steady state solutions, present traditional linear analysis of the problem 
and investigate analytically and numerically the influence of noise on the mean velocity of 
the front and pole dynamics. 

In Chapter 3 we obtain results of the linear stability analysis by the help of pole solutions. 

In Chapter 4 we use the result obtained for the channel case for analysis of the flame 
front propagation in the radial case 

In Chapter 5 we investigate asymptotic behavior of the poles in the complex plane for the 
Laplacian growth with the zero surface-tension in the case of periodic and no-flux boundary 
condition. 

Chapter 6 is a summary. 
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Chapter 2 



Pole-Dynamics in Unstable Front 
Propagation: the Case of the Channel 
Geometry 

2.1 Introduction 

The aim of this chapter is to examine the role of random fluctuations on the dynamics 
of growing wrinkled interfaces which are governed by non-linear equations of motion. We 
are interested in those examples for which the growth of a flat or smooth interface are 
inherently unstable. A famous example of such growth phenomena is provided by Laplacian 
growth patterns [1-3]. The experimental realization of such patterns is seen for example 
in Hele-Shaw cells [1] in which air or another low viscosity fluid is displacing oil or some 
other high viscosity fluid. Under normal conditions the advancing fronts do not remain 
flat; in channel geometries they form in time a stable finger whose width is determined 
by delicate effects that arise from the existence of surface tension. In radial geometry, the 
growth the interface forms a contorted and ramified fractal shape. A related phenomenon has 
been studied in a model equation for flame propagation which has the same linear stability 
properties as the Laplacian growth problem [9]. The physical problem in this case is that 
of pre-mixed flames which exist as self-sustaining fronts of exothermic chemical reactions in 
gaseous combustion. Experiments [10] on flame propagation in radial geometry show that 
the flame front accelerates as time goes on, and roughens with characteristic exponents. Both 
observations did not receive proper theoretical explanations. It is notable that the channel 
and radial growth are markedly different; the former leads to a single giant cusp in the moving 
front, whereas the latter exhibits infinitely many cusps that appear in a complex hierarchy 
as the flame front develops ( [11,19] and chapter 4). 

Analytic techniques to study such processes are available [38]. In the context of flame 
propagation [12,15,19,39], and in Laplacian growth in the zero surface-tension limit [6,35,36] 
one can examine solutions that are described in terms of poles in the complex plane. This 
description is very useful in providing a set of ordinary differential equations for the positions 
of the poles, from which one can deduce the geometry of the developing front in an extremely 
economical and efficient way. Unfortunately this description is not available in the case of 
Laplacian growth with surface tension, and this makes the flame propagation problem very 
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attractive. However, it suffers from one fundamental drawback. For the noiseless equation the 
pole-dynamics always conserves the number of poles that existed in the initial conditions. As 
a result there is a final degree of ramification that is afforded by every set of initial conditions 
even in the radial geometry, and it is not obvious how to describe the continuing self-similar 
growth that is seen in experimental conditions or numerical simulations. Furthermore, as 
mentioned before, at least in the case of flame propagation one observes [10] an acceleration 
of the flame front with time. Such a phenomenon is impossible when the number of poles is 
conserved. It is therefore tempting to conjecture that noise may have an important role in 
affecting the actual growth phenomena that are observed in such systems. In fact, the effect 
of noise on unstable front dynamics has not been adequately addressed in the literature. 
From the point of view of analytic techniques noise can certainly generate new poles even if 
the initial conditions had a flnite number of poles. The subject of pole dynamics with the 
existence of random noise, and the interaction between random fluctuations and deterministic 
front propagation are the main issues of this chapter. 

We opt to study the example of flame propagation rather than Laplacian growth, simply 
because the former has an analytic description in terms of poles also in the experimentally 
relevant case of finite viscosity. We choose to begin the study with channel geometry. The 
reason is that in radial geometry it is more difficult to disentangle the effects of external noise 
from those of initial conditions. After all, initially the system can contain infinitely many 
poles,very far away near infinity in the complex plane (and therefore having an infinitely 
small contribution to the interface). Since the growth of the radius changes the stability 
of the system, more and more of these poles might fall down to the real axis and become 
observable. In channel geometry the analysis of the effect of initial conditions is relatively 
straightforward, and one can understand it before focusing on the (more interesting) effects 
of external noise [12]. The basic reason for this is that in this geometry the noiseless steady 
state solution for the developed front is known analytically. As described in Section II, in a 
channel of width L the steady-state solution is given in terms of N{L) poles that are organized 
on a line parallel to the imaginary axis. It can be shown that for any number of poles in 
the initial conditions this is the only attractor of the pole dynamics. After the establishment 
of this steady state we can begin to systematically examine the effects of external noise on 
this solution. As stated before, in radial conditions there is no stable steady state with a 
finite number of poles, and the disentanglement of initial vs. external perturbations is less 
straightforward ( [19] and chapter 4). We show later that the insights provided in this chapter 
have relevance for radial growth as well as will be discussed in the sequel. 

We have a number of goals in this chapter. Firstly, after introducing the pole decom- 
position, the pole dynamics, and the basic steady state, we will present stability analysis of 
the solutions of the flame propagation problem in a channel geometry. It will be shown that 
the giant cusp solution is linearly stable, but non-linearly unstable. These results, which are 
described in Section III, can be obtained either by linearizing the dynamics around the giant 
cusp solutions in order to study the stability eigenvalues, or by examining perturbations in 
the form of poles in the complex plane. The main result of Section III is that there exists 
one Goldstone mode and two modes whose eigenvalues hit the real axis periodically when 
the system size L increases. Thus the system is marginally stable at particular values of 
L, and it is always nonlinearly unstable, allowing flnite size perturbations to introduce new 
poles into the system. This insight allows us to understand the relation between the system 
size and the effects of noise. In Section IV we discuss the relaxation dynamics that ensues 
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after starting the system with "small" initial data. We study the coarsening process that 
leads in time to the final solution of the giant cusp, and understand from this what are the 
typical time scales that exist in our dynamics. We offer in this Section some results of nu- 
merical simulations that are interpreted in the later sections. In Section V we focus on the 
phenomenon of acceleration of the flame front and its relation to the existence of noise. In 
noiseless conditions the velocity of the flame front in a finite channel is bounded [12]. This 
can be shown either by using the pole dynamics or directly from the equation of motion. We 
will present the results of numerical simulations where the noise is controlled, and show how 
the velocity of the flame front is affected by the level of the noise and the system size. The 
main results are: (i) Noise is responsible for introducing new poles to the system; (ii) For low 
levels of noise the velocity of the flame front scales with the system size with a characteristic 
exponent; (iii) There is a phase transition at a sharp (but system-size dependent) value of 
the noise-level, after which the behavior of the system changes qualitatively; (iv) After the 
phase transition the velocity of the flame front changes very rapidly with the noise level. In 
the last Section we remark on the implications of these observations for the scaling behavior 
of the radial growth problem, and present a summary and conclusions. 



2.2 Equations of Motion and Pole-decomposition in the 
Channel Geometry 



It is known that planar flames freely propagating through initially motionless homogeneous 
combustible mixtures are intrinsically unstable. It was reported that such flames develop 
characteristic structures which include cusps, and that under usual experimental conditions 
the flame front accelerates as time goes on. A model in 1 -|- 1 dimensions that pertains to 
the propagation of flame fronts in channels of width L was proposed in [9]. It is written in 
terms of position h{x, t) of the flame front above the x-axis. After appropriate rescalings it 
takes the form: 



The domain is < x < L, z/ is a parameter and we use periodic boundary conditions. The 
functional I[h{x,t)\ is the Hilbert transform which is conveniently deflned in terms of the 
spatial Fourier transform 



For the purpose of introducing the pole-decomposition it is convenient to rescale the domain 
to < ^ < 27r. Performing this rescahng and denoting the resulting quantities with the same 
notation we have 




(2.1) 




(2.2) 



I[h{k,t)] = \k\h{k,t) 



(2.3) 



dt 2L2 [ de 

+\l{h{6,t)} + l. 



dh{9,t) _ 1 \dh{e,t) 



+ 



V d'^h{e,t) 

L2 



(2.4) 
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In this equation L — L/27r. Next we change variables to u{9, t) = dh{9, t)/d9. We find 



(2.5) 



It is well known that the flat front solution of this equation is linearly unstable. The linear 
spectrum in /^-representation is 



ujk = \k\/L - vk'^/L^ . 
There exists a typical scale k^ax which is the last unstable mode 

k 

"'max 



(2.6) 



(2.7) 



Nonlinear effects stabilize a new steady-state which is discussed next. 

The outstanding feature of the solutions of this equation is the appearance of cusp-like 
structures in the developing fronts. Therefore a representation in terms of Fourier modes is 
very inefficient. Rather, it appears very worthwhile to represent such solutions in terms of 
sums of functions of poles in the complex plane. It will be shown below that the position of 
the cusp along the front is determined by the real coordinate of the pole, whereas the height 
of the cusp is in correspondence with the imaginary coordinate. Moreover, it will be seen 
that the dynamics of the developing front can be usefully described in terms of the dynamics 
of the poles. Following [12,19,38,39] we expand the solutions u{6,t) in functions that depend 
on N poles whose position Zj{t) = Xj{t) + iyjit) in the complex plane is time dependent: 



N 

u{9, t) = cot 

N 



9-zj{t) 



+ c.c. 



2sm[9 - Xj{t)] 



\ cosh[yj(t)] — cos[^ — Xj{t)\ 



(2.8) 



h{9,t) = 2z/^ln [cosh(?/j(t)) 



cos 



x,{t)) +C(t). 



(2.9) 



In (2.9) C(t) is a function of time. The function (2.9) is a superposition of quasi-cusps 
(i.e. cusps that are rounded at the tip). The real part of the pole position (i.e. Xj) is the 
coordinate (in the domain [0, 27r]) of the maximum of the quasi-cusp, and the imaginary part 
of the pole position (i.e |/j) is related to the depth of the quasi-cusp. As yj decreases the 
depth of the cusp increases. As — > the depth diverges to infinity. Conversely, when 
yj — > oo the depth decreases to zero. 

The main advantage of this representation is that the propagation and wrinkling of the 
front can be described via the dynamics of the poles. Substituting (2.8) in (2.5) we derive 
the following ordinary differential equations for the positions of the poles: 



dz 
dt 



2N 

V ^ cot 



(2.10) 
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We note that in (2.8), due to the complex conjugation, we have 2N poles which are arranged 
in pairs such that for j < N Zj^j\f = Zj. In the second sum in (2.8) each pair of poles 
contributed one term. In Eq.(2.10) we again employ 2N poles since all of them interact. We 
can write the pole dynamics in terms of the real and imaginary parts Xj and yj. Because of 
the arrangement in pairs it is sufficient to write the equation for either i/j > or for i/j < 0. 
We opt for the first. The equations for the positions of the poles read 



dx- ^ 
df ^ 



sm[x-i 



[cosh{yj - yk) 



(2.11) 



cos{xj - Xk)] ^ + [cosh{yj + yk) - cos{xj - Xk)] ^ 
>dyj ^ f smh{yj 



dt 



Vk) 



k=i^k^j ^cosh(|/j- - yk) - cos{xj - Xk) 
sinh(?/j + yk) 



cosh(?/j + yk) - cos{xj - Xk] 



+ ucoth.{yj) — L. 



(2.12) 



We note that if the initial conditions of the differential equation (2.5) are expandable in a 
finite number of poles, these equations of motion preserve this number as a function of time. 
On the other hand, this may be an unstable situation for the partial differential equation, 
and noise can change the number of poles. This issue will be examined at length in Section 
2.5. 



2.3 Linear Stability Analysis in Channel Geometry 

In this section we discuss the linear stability of the TFH-cusp solution. To this aim we first 
use Eq.(2.8) to write the steady solution Us{0) in the form: 

m\ 2sin[^-a;,] 

'U's{(')=^y — n — r-^ FTT (2-13) 

' fr[ cosh[yj] - cos[^ - X,] ' ^ ' 

where Xs is the real (common) position of the stationary poles and yj their stationary imagi- 
nary position. To study the stability of this solution we need to determine the actual positions 
yj. This is done numerically by integrating the equations of motion for the poles starting 
from poles in initial positions and waiting for relaxation. Next one perturbs this solution 
with a small perturbation 0(6*, t): u{9,t) = Us{6) + (l){9,t) . Linearizing the dynamics for 
small (f) results in the equation of motion 

= ^ ldo[us(e)m t)] + ud'.m t) 

+ Tl{mt)) . (2.14) 
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Figure 2.1: The first 10 highest eigenvalues of the stabihty matrix with v = n/5, multiphed 
by the square of the system size vs. the system size L. Note that all the eigennvalues 
oscillate around fixed values in this presentation, and that the highest two eigenvalues hit 
zero periodically. 

2.3.1 Fourier decomposition and eigenvalues 

The hnear equation can be decomposed in Fourier modes according to 

oo 

<l>{e,t) = Yl Mt)e"'' (2.15) 

fe=— oo 

oo N 

Us{e) = -2ui Esign(fc)e~l'=l^^e^'=^ (2.16) 

fc=— oo j=l 

In these sums the discrete k values run over all the integers. Substituting in (2.14) we get 
the equations: 

#fe(^)) ? /A /o 1'7^ 

dt 

where Ukn is a infinite matrix whose entries are given by 

akk = ^-X^^' (2-18) 

k ^ 
akn = — sign(fc-n)(2z/^e-l'=-'^l^0 k^n. (2.19) 

To solve for the eigenvalues of this matrix we need to truncate it at some cutoff A;- vector k*. 
The choice of k* can be based on the linear stabihty analysis of the flat front. The scale kjnax: 
cf. (2.7), is the largest k which is still linearly unstable. We must choose k* > k^ax ^-nd test 
the choice by the converegence of the eigenvalues. The chosen value of k* in our numerics 
was 4:kmax- The results for the low order eigenvalues of the matrix akn that were obtained 
from a converged numerical calculation are presented in Fig. 2.1. 

The eigenvalues are multiphed by and are plotted as a function of L. We order the 
eigenvalues in decreasing order and denote them as | Aq |<| Ai |<| A2 | 



14 



Fig 2.1 contains a strange result on the positive eigenvalues at large L. One of methods 
to check some numerical result is to do analytic investigation. For example, in Chapter 
3 we make detailed analytic investigation for the numerical result on Fig. 2.1 and obtain 
that all eigenvalues are not positive. Indeed, two types of modes exists. The first one is 
connected to the displacement of poles in the giant cusp. Because of the pole attraction 
the giant cusp is stable with respect to the longitudinal displacement of poles and so the 
correspondent eigenvalues are not positive. For the transversal displacement the Lyapunov 
function exists and so the giant cusp is stable with respect to the transversal displacement 
and the correspondent eigenvalues are not positive. The second type of modes is connected to 
additional poles. These poles go to infinity because of the repulsion from the giant cusp poles 
N{L). So the correspondent eigenvalues are also not positive. So the positive eigenvalues at 
large L are a numerical artifact. 

The figure ofi^crs a number of qualitative observations: 

1. There exists an obvious Goldstone or translational mode u'g{0) with eigenvalue Aq = 0, 
which is shown with rhombes in Fig.2.1. This eigenmode stems from the Galilean 
invariance of the equation of motion. 

2. The eigenvalues oscillate periodically between values that are L-independent in this 
presentation (in which we multiply by L^). In other words, up to the oscillatory behavior 
the eigenvalues depend on L like L~'^. 

3. The eigenvalues Ai and A2, which are represented by squares and circles in Fig.2.1, 
hit zero periodically. The functional dependence in this presentation appears almost 
piecewise linear. 

4. The higher eigenvalues also exhibit similar qualitative behaviour, but without reaching 
zero. We note that the solution becomes marginally stable for every value of L for which 
the eigenvalues Ai and A2 hit zero. The L^^ dependence of the spectrum indicates that 
the solution becomes more and more sensitive to noise as L increases. 



2.3.2 Qualitative understanding using pole-analysis 

The most interesting qualitative aspects are those enumerated above as item 2 and 3. To 
understand them it is useful to return to the pole description, and to focus on Eq.(l.ll). 
This equation describes the dynamics of a single far-away pole. We remarked before that 
this equation shows that for fixed L the stable number of poles is the integer part (1.10). 
Define now the number a, < a < 1, according to 



a — 



1(1-1). (2.20) 



Using this number we rewrite Eq.(l.ll) as 

As L increases, a oscillates piecewise linearly and periodically between zero and unity. This 
shows that a distant pole which is added to the giant cusp solution is usually repelled to 
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Figure 2.2: Comparison of the numerically determined highest 4 eigenvalues of the stability 
matrix with the prediction of the pole analysis. The eigenvalues of the stabilit y matrix are 
:Ao, Ai, A2 and A3. The pole analysis (solid line) provides a qualitative understanding of the 
stability, and appears to overlap with the high est eigenvector over half of the range, and 
with the fourth eigenvalue over the other half. 

infinity except when a hits zero and the system becomes marginally unstable to the addition 
of a new pole. 

To connect this to the linear stability analysis we note from Eq.(2.8) that a single far-away 
pole solution (i.e with y very large) can be written as 

u{e, t) = 4z/e-^(*) sin(^ - x{t)) . (2.22) 

Suppose that we add to our giant cusp solution a perturbation of this functional form. 
^Prom Eq.(2.21) we know that y grows linearly in time, and therefore this solution decays 
exponentially in time. The rate of decay is a linear eigenvalue of the stability problem, and 
from Eq.(2.21) we understand both the dependence and the periodic marginality. We 
should note that this way of thinking gives us a significant part of the L dependence of the 
eigenvalues, but not all. The variable a is rising from zero to unity periodically, but after 
reaching unity it hits zero instantly. Accordingly, if the highest non zero eigenvalue were 
fully determined by the pole analysis, we would expect this eigenvalue to behave as the solid 
line shown in Fig. 2. 2. 

The actual highest eigenvalue computed from the stability matrix is shown in rhombes 
connected by dotted line. It is clear that the pole analysis gives us a great deal of qualitative 
and quantitative understanding, but not all the features agree. 

2.3.3 Dynamics near marginality 

The discovery of marginality at isolated values of L poses questions regarding the fate of 
poles that arc added at very large y's at certain ,x-positions. We will argue now that when 
the system becomes marginally stable, a new pole can be added to those existing in the giant 
cusp. We remember that these poles have a common 9 position that we denote as 6 — 9c. 
The fate of a new pole added at infinity depends on its 9 position. If the position of the new 
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pole is again denoted as Ua, and oo ^ i/a ^ Umax, can see from Eq.(2.12) that di/a/dt is 
maximal when 6a = 0^, whereasit is minimal when 6a — = vr. This follows from the fact 
that the cosine term has a value +1 when 9a = 9c and a value —1 when 9a — 9c = tt. For large 
y differences the terms in the sum take on their minimal value when the cos term is —1 and 
their maximal values at +1. For infinitely large i/a the equation of motion is (1-11) which is 
independent of 9a- Since the RHS of this equation becomes zero at marginality, we conclude 
that for very large but finite ya dya/ dt changes sign from positive to negative when 9a — 9^ 
changes from zero to tt. The meaning of this observation is that the most unstable points in 
the system are those points which are furthest away from the giant cusp. It is interesting to 
discuss the fate of a pole that is added to the system at such a position. From the point of 
view of the pole dynamics = + is an unstable fixed point for the motion along the 9 
axis. The attraction to the giant cusp exactly vanishes at this point. If we start with a pole 
at a very large ya close to this value of 9 the down-fall along the y coordinate will be faster 
than the lateral motion towards the giant cusp. We expect to see therefore the creation of 
a small cusp at 9 values close to tt that precedes a later stage of motion in which the small 
cusp moves to merge with the giant cusp. Upon the approach of the new pole to the giant 
cusp all the existing poles will move up and the furthest pole at ymax will be kicked off to 
infinity. We will later explain that this type of dynamics occurs in stable systems that are 
driven by noise. The noise generates far away poles (in the imaginary direction) that get 
attracted around ^ = + tt to create small cusps that run continuously towards the giant 
cusp. 

2.3.4 Excitable System. 

The intuition gained so far can be used to discuss the issue of stability of a stable system 
to larger perturbations. In other words, we may want to add to the system poles at finite 
values of y and ask about their fate. We first show in this subsection that poles whose initial 
y value is below ymax ~ ^og{L'^ /v^) will be attracted towards the real axis. The scenario is 
similar to the one described in the last paragraph. 

Suppose that we generate a stable system with a giant cusp at = with poles distributed 
along the y axis up to ymax- We know that the sum of all the forces that act on the upper 
pole is zero. Consider then an additional pole inserted in the position (tt, ymax)- It is obvious 
from Eq.(2.12) that the forces acting on this pole will pull it downward. On the other hand 
if its initial position is much above ymax the force on it will be repulsive towards infinity. We 
see that this simple argument identifies ymax as the typical scale for nonlinear instability. 

Next we estimate ymax and interpret our result in terms of the amplitude of a perturbation 
of the flame front. We explained that uppermost pole's position fluctuates between a minimal 
value and inflnity as L is changing. We want to estimate the characteristic scale of the 
minimal value of ymax{L). To this aim we employ the result of rcf. [12] regarding the stable 
distribution of pole positions in a stable large system. The parametrization of [12] differs 
from ours; to go from our parametrization in Eq.(2.5) to theirs we need to rescale u by 
and t by L. The parameter u in their parameterization is v/L in ours. According to [12] the 
number of poles between y and y + dy is given by the p{y)dy where the density p{y) is 

= Ain[coth(|t/|/4)] . (2.23) 
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To estimate the minimal value of Umax we require that the tail of the distribution p[y) 
integrated between this value and infinity will allow one single pole. In other words, 



/ dypiy) « 1 . (2.24) 

Umax 



Expanding (2.23) for large y and integrating explicitly the result in (2.24) we end up with 
the estimate 

Vmax f« 2 In — (2.25) 

For large L this result is ymax ~ In(^). If we now add an additional pole in the position 
(S^ymax) this is equivalent to perturbing the solution u{9,t) with a function i/e"^'""'^ sin(^), 
as can be seen directly from (2.8). We thus conclude that the system is unstable to a 
perturbation larger than 

u{e) sm{e) / L"^ . (2.26) 

This indicates a very strong size dependence of the sensitivity of the giant cusp solution 
to external perturbations. This will be an important ingredient in our discussion of noisy 
systems. 



2.4 Initial Conditions, Pole Decomposition and Coars- 
ening 

In this section we show first that any initial conditions can be approximated by pole de- 
composition. Later, we show that the dynamics of sufficiently smooth initial data can be 
well understood from the pole decomposition. Finally we employ this picture to describe 
the inverse cascade of cusps into the giant cusp which is the final steady state. By inverse 
cascade we mean a nonlinear coarsening process in which the small scales coalesce in favor 
of larger scales and finally the system staturates at the largest available scale [40] . 



2.4.1 Pole Expansion: General Comments 

The fundamental question is how many poles are needed to describe any given initial con- 
dition. The answer, of course, depends on how smooth are the initial conditions. Suppose 
also that we have an initial function u{9, t — Q) that is 27r-periodic and which at time t — Q 
admits a Fourier representation 

oo 

u{e) = Y,AkSin{ke + <Pk) , (2.27) 
fe=i 

with Ak > Q for all k. Suppose that we want to find a pole-decomposition representation 
Up{9) such that 

\up{e)-u{e)\ < e for every 9 , (2.28) 

where e is a given wanted accuracy. If u{9) is differentiable we can cut the Fourier expansion 
at some finite k — K knowing that the remainder is smaller than, say, e/2. Choose now a 
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large number M and a small number A <^ 1/M and write the pole representation for Up{9) 
as 



fc=l p=0 



cosh [k{yk + pA)] — cos{k6 + (pk) 



To see that this representation is a particular form of the general formula (2.8) We use the 
following two identities 

f^e-^'smxk = - /^"^ , (2.30) 

2 cosh i — cos X 

■ ( . ,K-l Ky y 

} ^ sm(x + ky) — sm(x H — y) sm -— cosec- . {^■•^^) 

k=o 2 2 2 



Prom these follows a third identity 



^1 2sin(a;-^ + 



Q cosh y — cos {x — ^ + (f)) 

2Ksm{Kx + 0) 
cosh Ky — cos {Kx + (f)) 



(2.32) 



Note that the LHS of (2.32) is of the form (2.8) with K poles whose positions are all on the 
line yj = y and whose Xj are on the lattice points 27rj/ K — (f). On the other hand every term 
in (2.29) is of this form. 

Next we use (2.30) to rewrite (2.29) in the form 

K M-l oo 

M^) = E E E 4A;e-"*^(^'=+P^) sin {nk9 + ncj)k) • (2.33) 

k=l p=0 n=l 

Exchanging order of summation between n and p we can perform the geometric sum on p. 
Denoting 

M-l 1 _ -MfcnA 

bn,k - E e-"''' - 1 _ ,-knA ' (2-34) 



we find 



p=0 



K oo 

= EE4fcVfee-"*^^'sin(nA;e + n(/)fc) 

fc=ln=l 
K oo 

= E E 4A:6„,fee-"'=2''= sin {nkO + n^^) 

A:=ln=2 

+ E ^khx^ke'^y^ sin (A;^ + 0^) . (2.35) 
fc=i 



Compare now the second term on the RHS of (2.35) with (2.27). We can identify 



(2.36) 
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The first term can be then bound from above as 



K CO 

J2 E 4A;6„,fee-"^^'= sin {nk9 + ncp^) 

k=l n=2 

K CO 



(2.37) 



k=l n=2 



A, 



sin [nkO + n<pk) 



The sine function and the factor (4i^)^ can be replaced by unity and we can bound the 
RHS of (2.37) by 



K oo 

EE 

k=l n=2 



A, 



K 



bn,k < E E 

k=l n=l 



A, 



(2.38) 



where we have used the fact that bn^k < which follows directly from (2.34). Using now 
the facts that bi^K < for every k < K and that Aj^ is bounded by some finite C since it 
is a Fourier coefficient, we can bound (2.38) by C'^K/{bi^K — C). Since we can select the free 
parameters A and M to make bi K as large as we want, we can make the remainder series 
smaller in absolute value than e/2. 

The conclusion of this demonstration is that any initial condition that can be represented 
in Fourier series can be approximated to a desired accuracy by pole-decomposition. The 
number of needed poles is of the order x M. Of course, the number of poles thus 
generated by the initial conditions may exceed the number N{L) found in Eq.(l.lO). In such 
a case the excess poles will move to infinity and will become irrelevant for the short time 
dynamics. Thus a smaller number of poles may be needed to describe the state at larger 
times than at i = 0. We need to stress at this point that the pole decomposition is over 
complete; for example, if there is exactly one pole at t = and we use the above technique 
to reach a pole decomposition we would get a large number of poles in our representation. 



2.4.2 The initial stages of the front evolution: the exponential 
stage and the inverse cascade 

In this section wc employ the connection between Fourier expansion and pole decomposition 
to understand the initial exponential stage of the evolution of the flame front with small 
initial data u{6,t = 0). Next we employ our knowledge of the pole interactions to explain 
the slow dynamics of coarsening into the steady state solution. 

Suppose that initially the expansion (2.27) is available with all the coefficients A^ <^ 1. 
We know from the linear instability of the flat flame front that each Fourier component 
changes exponentially in time according to the linear spectrum (2.6). The components with 
wave vector larger than (2.7) decrease, whereas those with lower wave vectors increase. The 
fastest growing mode is kc — L/2y. In the linear stage of growth this mode will dominate 
the shape of the flame front, i.e. 

u{e,t) ^ Ak^e'^'^-' sm{kc9) . (2.39) 

Using Eq.(2.32) for a large value of y (which is equivalent to small A^J we see that to the order 
of 0(A|, ) (2.39) can be represented as a sum over L/2u poles arranged periodically along 
the 6 axis. Other unstable modes will contribute similar arrays of poles but at much higher 
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values of y, since their amplitude is exponentially smaller. In addition we have nonlinear 
corrections to the identification of the modes in terms of poles. These corrections can be 
again expanded in terms of Fourier modes, and again identified with poles, which will be 
further away along the y axis, and with higher frequencies. To see this one can use Eq.(2.35), 
subtract from Up{6) the leading pole representations, and reexpand in Fourier series. Then 
we identify the leading order with double the number of poles that are situated twice further 
away along the y axis. 

We note that even when all the unstable modes are present, the number of poles in the first 
order identification is finite for finite L, since there are only L/i/ unstable modes. Counting 

the number of poles that each mode introduces we get a total number of {l/^^ poles. 
The number L/2v oi poles which are associated with the most unstable mode is precisely the 
number allowed in the stable stationary solution, cf.(l.lO). When the poles approach the real 
axis and cusps begin to develop, the linear analysis no longer holds, but the pole description 
does. 

We now describe the qualitative scenario for the establishment of the steady state. Firstly, 
we understand that all the poles that belong to less unstable modes will be pushed towards 
infinity. To see this think of the system at this stage as an array of uncoupled systems with 
a scale of the order of unity. Each such system will have a characteristic value of y. As 
we discussed before poles that are further away along the y axis will be pushed to infinity. 
Therefore the system will remain with the L/2v poles of the most unstable mode. The net 
effect of the poles belonging to the (nonlinearly) stable modes is to destroy the otherwise 
perfect periodicity of the poles of the unstable mode. The see the effect of the higher order 
correction to the pole identification we again recall that they can be represented as further 
away poles with higher frequencies, whose dynamics is similar to the less unstable modes 
that were just discussed. They do not become more relevant when time goes on. 

Once the poles of the stable modes get sufficiently far from the real axis, the dynamics of 
the remaining poles will begin to develop according to the interactions that are directed along 
the real axis. These interactions are much weaker and the resulting dynamics occur on much 
longer time scales. The qualitative picture is of an inverse cascade of merging the 6 positions 
of the poles. We note that the system has a set of unstable fixed points which are 'cellular 
solutions' described by a periodic arrangement of poles along the real axis with a frequency k. 
These fixed points are not stable and they collapse, under perturbations, with a characteristic 
time scale (that depends on k) to the next unstable fixed point at k' = k/2. This process 
then goes on indefinitely until A; ~ 1/L i.e. we reach the giant cusp, the steady-state stable 
solution [40]. 

This scenario is seen very clearly in the numerical simulations. In Fig.2.3 we show 
the time evolution of the flame front starting from small white-noise initial conditions. 
The bottom curve pertains to the earliest time in this picture, just after the fast exponential 
growth, and one sees clearly the periodic array of cusps that form. The successive images 
show the progress of the fiame front in time, and one observes the development of larger scales 
with deeper cusps that represent the partial coalescence of poles onto the same 9 positions. 
In Fig.2.4 

we show the width and the velocity of this front as a function of time. One recognizes the 
exponential stage of growth in which the L/2i/ poles approach the 9 axis, and then a clear 
cross-over to much slower dynamics in which the effective scale in the system grows with a 
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Figure 2.3: The inverse cascade process of coarsening that occurs after preparing the system 
with random, small initial conditions. One sees that at successive times the typical scale 
increases until the giant cusp forms, and attracts all the other side-poles. The effect of the 
existing numerical additive noise is to introduce poles that appear as side cusps that are 
continuously attracted to the giant cusp. This effect is obvious to the eye only after the 
typical scale is sufficiently large, as is seen in the last time (see text for further details). 




Figure 2.4: log-log plots of the front velocity (lower curve) and width (upper curve) as a 
function of time in the inverse cascade process seen in Fig. 2. 3 in a system o f size 2000 and 
z/ = 1. Both quantities exhibit an initial exponential growth that turns to a power law growth 
(after t ^ 30). The velocity is constant after this time, and the width increases like t^. Note 
that at the earliest time there is a slight decrease in the velocity; this is due to the decay of 
linearly stable modes that exist in random initial conditions. 
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slower rate. 

The slow dynamics stage can be understood qualitatively using the previous interpretation 
of the cascade as follows: if the initial number of poles belonging to the unstable mode is 
L/2v, the initial effective linear scale is 2i/. Thus the first step of the inverse cascade will be 
completed in a time scale of the order of 2v. At this point the effective linear scale doubles 
to 4z/, and the second step will be completed after such a time scale. We want to know what 
is the typical length scale It seen in the system at time t. The definition of front width is 

It — \j i /o^[/i(a;, t) — h]'^dx,h = i Jq h{x, i)dx. The typical width of the system at this stage 
will be proportional to this scale. 

Denote the number of cascade steps that took place until this scale is achieved by Si. The 
total time elapsed, t{lt) is the sum 

t{lt)r^Y.T. (2.40) 

i=l 

The geometric sum is dominated by the largest term and we therefore estimate t{lt) ~ U- 
We conclude that the scale and the width are linear in the time elapsed from the initial 
conditions (It ~ t'', ^ = 1). In noiseless simulations we find (see Fig.2.4) a value of ^ which 
is C ~ 0.95 ±0.1. 

2.4.3 Inverse cascade in the presence of noise 

An interesting consequence of the discussion in the last section is that the inverse cascade 
process is an effective "clock" that measures the typical time scales in this system. For future 
purposes we need to know the typical time scales when the dynamics is perturbed by random 
noise. To this aim we ran simulations following the inverse cascade in the presence of external 
noise. The main result that will be used in later arguments is that now the appearance of a 
typical scale It occurs not after time t, but rather according to 

k^t^ , C~ 1-2 ±0.1. (2.41) 

The numerical confirmation of this law is exhibited in Fig. 2. 5 . 

We also find that the front velocity in this case increases with time according to 

~ r , 7 0.48 ± 0.05 . (2.42) 

This result will be related to the acceleration of the fiame front in noisy simulations, as will 
be seen in the next Sections. 

2.5 Acceleration of the Flame Front, Pole Dynamics 
and Noise 

A major motivation of this Section is the observation that in radial geometry the same 
equation of motion shows an acceleration of the flame front. The aim of this section is to 
argue that this phenomenon is caused by the noisy generation of new poles. Moreover, it is 
our contention that a great deal can be learned about the acceleration in radial geometry by 
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Figure 2.5: The same as Fig. 2. 4 but with additive random noise for a system of size 1000, 
u = 0.1 and / = 10~^^. The velocity does not saturate now, and the exponent ( characterizing 
the increase of the width with time changes to ( = 1.2±0.1. The velocity increases in time 
like P with 7 ^ 0.48 ± 0.04. 

considering the effect of noise in channel growth. In Ref. [12] it was shown that any initial 
condition which is represented in poles goes to a unique stationary state which is the giant 
cusp which propagates with a constant velocity f = 1/2 up to small 1/L corrections. In light 
of our discussion of the last section we expect that any smooth enough initial condition will 
go to the same stationary state. Thus if there is no noise in the dynamics of a finite channel, 
no acceleration of the flame front is possible. What happens if we add noise to the system? 

For concreteness we introduce an additive white-noise term 77(6', t)to the equation of mo- 
tion (2.5) where 

r]{e,t)=Y.r]k{t)exp{ike) , (2.43) 

k 

and the Fourier amplitudes rjk are correlated according to 

< VkitWAt') >= ^h,k'6it - t') . (2.44) 

We will first examine the result of numerical simulations of noise-driven dynamics, and later 
return to the theoretical analysis. 

2.5.1 Noisy Simulations 

Previous numerical investigations [11,13] did not introduce noise in a controlled fashion. We 
will argue later that some of the phenomena encountered in these simulations can be ascribed 
to the (uncontrolled) numerical noise. We performed numerical simulations of Eq.(2.5 using 
a pseudo-spectral method. The time-stepping scheme was chosen as Adams-Bashforth with 
2nd order presicion in time. The additive white noise was generated in Fourier-space by 
choosing r/^ for every k from a flat distribution in the interval [— y^2^, •\/2^]- We examined 
the average steady state velocity of the front as a function of L for fixed / and as a function 
of / for fixed L. We found the interesting phenomena that are summarized here: 

1. In Fig.2.7 we can see two different regimes of the behavior of the average velocity v as 
a function of the noise /"'^ for the fixed system size L. For the noise / smaller then 
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Figure 2.6: The dependence of the average velocity v on the system size L for /^'^ = 0, 2.7 x 
10-^ 2.7 X 10-^ 2.7 X 10-^ 2.7 X 10-3, 2.7 x IQ-^, 2.7 x 10-\ 0.5, 1.3, 2.7. 
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Figure 2.7: The dependence of the average velocity v on the noise /°'^ for L=10, 40, 80. 
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Figure 2.8: Typical flame fronts for / < fcr where the system is sufficientl y smaU not to be 
terribly affected by the noise. The effect of noise in this regime is to add additional small 
cusps to the giant cusp. In figures a-d we present fronts for growing system sizes L = 10, 20, 40 
and 80 respectively, u = 0.1 . One can observe that when the system size grows there are 
more cusps with a more complex structure. 

same fixed value fcr 

vr^ f . (2.45) 

For these values of / this dependence is very weak, and ^ 0.02. For the large values 
of / the dependence is much stronger 

2. In Fig. 2. 6 we can see the growth of the average velocity v as a, function of the system 
size L. After some values of L we can see saturation of the velocity. For regime f < fcr 
the growth of the velocity can be written as 

vr^L'', 0.35 ± 0.03 . (2.46) 

3. In Fig. 2. 8 and Fig.2.9 we can see fiame fronts for f < far and / > fcr- 

2.5.2 Calculation of the Number of Poles in the System 

The interesting problem that we would like to solve here to better understand the dynamics 
of poles, is to determine those that exist in our system outside the giant cusp. This can be 
done by calculating the number of cusps (points of minimum or infiexional points) and their 
position on the interval ^ : [0, 27r] in every moment of time and drawing the positions of the 
cusps like functions of time, see Fig. 2.10. In this picture we can see the x-positions of all 
cusps in the system as a function of time. 

We have assumed that our system is in a "quasi-stable" state most of the time, i.e. every 
new cusp that appears in the system includes only one pole. Using pictures obtained in this 
way we can find: 

1. The mean number of poles in the system. By calculating the number of cusps in some 
moment of time and by investigating the history of every cusp (except the giant cusp), 
i.e. how many initial cusps take part in formating this cusp, and after averaging the 
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Figure 2.9: A typical flame front for / > fcr- The system size is 160. This is sufficient to cause 
a qualitative change in the appearance of the flame front: the noise introduces signiflcant 
levels of small scales structure in addition to the cusps. 
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Figure 2.10: The dependence of the cusps positions on timea. L = 80i/ = 0.1/ = 9xl0 ^ 
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Figure 2.11: The dependence of the pole number in the unit time dN/dt on the noise /°'^. 
i/ = 0.1 L = 80 

number of poles found with respect to different moments of time, we can find the 
mean number of poles that exist in our system outside the giant cusp. Let us denote 
this number by 6N. There are four regimes that can be defined with respect to the 
dependence of this number on the noise /: 

(i) Regime I: Such httle noise that no new cusps exist in our system outside the giant 
cusp; 

(ii) Regime II: Strong dependence of the pole number 5N on the noise /; 

(iii) Regime III: Saturation of the pole number 5N on the noise /, so that this number 
depends very little on the noise (Fig. 2.12); 

p0.03 



5N ~ (2.47) 



The saturated value of SN is defined by next formula (Fig. 2.14, Fig. 2.16) 

SN ^ N{L)/2^^^ (2.48) 
where N(L) ^ is the number of poles in the giant cusp. 

(iv) Regime IV: We again see a strong dependence of the pole number SN on the noise 
/ (Fig. 2.12); 

5N - f -^ (2.49) 

Because of the numerical noise we can see in most of the simulations only regime III 
and IV. In the future if no new evidence is seen we will discuss regime III. 
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Figure 2.12: The dependence of the excess pole number 5N on the noise f . v = 0.1 
L = 40,80. 
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Figure 2.13: The dependence of the pole number in the unit time dN/dt on the system size 
L. z/ = 0.1 f-^ = 9 X 10"^ 
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Figure 2.14: The dependence of the excess pole number 5N on the system size L. u 
/0-5 = 9 X 10-6 



0.1 



10 



Figure 2.15: The dependence of the pole number in the unit time dN/dt on the parameter 
i/. L = 80 1/ = 0.1. 
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Figure 2.16: The dependence of the excess pole number 6N on the the parameter i>. L = 80 
i/ = 0.1. 

2. By calculating the new cusp number that appears in the system in the unit time we 
can find the number of poles that appear in the system in the unit time In regime 
III (Fig. 2.11) 



dN 0.03 

dt ■' 

The dependence on L and u is defined by (Fig. 2.13 and Fig. 2.15) 

dN 



(2.50) 



-0.8 



dt 



dN 



[2.51) 



dt 1/2 

In regime IV, the dependence on the noise is defined by the following: (Fig. 2.11) 

dN 



(2.52) 



dt 



f 



0.1 



(2.53) 



2.5.3 Theoretical Discussion of the Effect of Noise 



The Threshold of Instability to Added Noise. Transition from regime I to regime 
II 

First we present the theoretical arguments that explain the sensitivity of the giant cusp 
solution to the effect of added noise. This sensitivity increases dramatically with increasing 
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the system size L. To see this we use again the relationship between the hnear stabihty 
analysis and the pole dynamics. 

Our additive noise introduces perturbations with all /c-vectors. We showed previously 
that the most unstable mode is the k — 1 component Aism{9). Thus the most effective 
noisy perturbation is r)i sin(^) which can potentially lead to a growth of the most unstable 
mode. Whether or not this mode will grow depends on the amplitude of the noise. To 
see this clearly we return to the pole description. For small values of the amplitude Ai we 
represent Aism{6) as a single pole solution of the functional form i/e~^sin6'. The y position 
is determined from y — — log \Ai\/iy, and the ^-position is ^ = tt for positive Ai and ^ = for 
negative Ai. Prom the analysis of Section III we know that for very small Ai the fate of the 
pole is to be pushed to infinity, independently of its 6 position; the dynamics is symmetric 
in Ai — Ai when ij is large enough. On the other hand when the value of Ai increases the 
symmetry is broken and the 9 position and the sign of Ai become very important. U Ai > 
there is a threshold value of y below which the pole is attracted down. On the other hand 
if Ai < 0, and = the repulsion from the poles of the giant cusp grows with decreasing y. 
We thus understand that qualitatively speaking the dynamics of Ai is characterized by an 
asymmetric "potential" according to 

V{A^) = \Al-aAl + ... . (2.55) 

/,From the linear stability analysis we know that A ~ J^/L'^, cf. Eq.(l.ll). We know further 
that the threshold for nonlinear instability is at Ai pa i/^/L"^, cf. Eq(2.26). This determines 
that value of the coefficient a fa 2/3i/^. The magnitude of the "potential" at the maximum is 

V{A^a.) ~ i^'/L^ . (2.56) 

The effect of the noise on the development of the mode Ai sin 9 can be understood from the 
following stochastic equation 

i. = -^+mW. (2.57) 

It is well known [41] that for such dynamics the rate of escape R over the "potential" barrier 
for small noise is proportional to 

i?~ -^exp-'^'/^^' . (2.58) 

The conclusion is that any arbitrarily tiny noise becomes effective when the system size 
increase and when v decreases. If we drive the system with noise of amplitude ^ the system 
can always be sensitive to this noise when its size exceeds a critical value Lc that is determined 
by f /Lc ~ u'^^/L^. This formula defines transition from regime I (no new cusps) to regime 
II. For L > Lf. the noise will introduce new poles into the system. Even numerical noise in 
simulations involving large size systems may have a macroscopic infiuence. 

The appearance of new poles must increase the velocity of the front. The velocity is 
proportional to the mean of {u/ Vf. New poles distort the giant cusp by additional smaller 
cusps on the wings of the giant cusp, increasing u^. Upon increasing the noise amplitude 
more and more smaller cusps appear in the front, and inevitably the velocity increases. This 
phenomenon is discussed quantitatively in Section 2.5. 



32 



7.0 



6.5 



CO 



^^6.0 



5.5 



0.0 



50.0 . 100.0 



150.0 



Figure 2.17: The dependence of the normahzed amphtude A^axL /i^ on the system size L 



Numerical verifing of the asymmetric "potential" form and dependence of the 
noise on Lc 

Prom the equations of the motion for poles we can find the distribution of poles in the giant 
cusp [12]. If we know the distribution of poles in the giant cusp we can then find the form 
of the "potential" and verify numerically expressions for values A, A^ax and ^g^^-* discussed 
previously. The connection between amplitude Ai and the position of the pole y is defined 
by Ai = 4z/e^^ and the connection between the potential function '^^^^^ and the position 

of the pole y is defined by formula ^g^^-* = 4z/^e^^, where ^ can be determined from the 
equation of the motion of the poles. We can find A^ax as the zero-point of ^g^^-* and A 
can be found as ^^-^tt^ for Ai = 0. Numerical measurements were made for the set of 
values L = 2ni', where n is a integer and n > 2. For our numerical measurements we use the 
constant u — 0.005 and the variable L, where L changes in the interval [1,150], or variable u 
that changes in the interval [0.005,0.05] and the constant L — 1. The results obtained follow: 

1. ^"'°3 as a function of L is almost a constant. (Fig. 2.17) 

2. ^"""i^^ as a function of u is almost a constant. (Fig. 2.18) 

3- ^^7" as a function of L is almost a constant. {An(^l) is defined by the position of the 
upper pole.) (Fig. 2.19) 

4. f"'""' as a function of v is almost a constant. (Fig. 2.20) 

5. The value of ^ as a function of L is a constant (Fig. 2.21). 

6. The value of - — as a function of is a constant ( Fig. 2.22 ). 
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Figure 2.18: The dependence of the normahzed amphtude AmaxL'^/i^^ on the parameter u. 
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Figure 2.19: The relationship between the amphtude defined by the minimum of the potential 
Amax and the amplitude defined by the position of the upper pole Ajv(l) as a function of the 
system size L. 
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Figure 2.20: The relationship between the amphtuda defined by the minimum of the potential 
Aynax and the amplitude defined by the position of the upper pole Ajv(l) as a function of the 
parameter i/. 
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Figure 2.21: The dependence of the normalized parameter AL^/i/ on the system size L. 
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Figure 2.22: The dependence of the normaUzed parameter XL^/u on the parameter u. 




We also verify the boundary between regime I (no new cusps) and regime II (new cusps 
appear). Fig. 2.23 shows the dependence of on L^. We can see that //L^ ~ ^/L^c- These 
results are in good agreement with the theory. 

The Noisy Steady State and its Collapse with Large Noise and System Size 

In this subsection we discuss the response of the giant cusp solution to noise levels that are 
able to introduce a large number of excess poles in addition to those existing in the giant 
cusp. We will denote the excess number of poles by 5N . The first question that we address 
is how difficult is it to insert yet an additional pole when there is already a given excess 
8N . To this aim we estimate the effective potential Vsn{^i) which is similar to (2.55) but is 
taking into account the existence of an excess number of poles. A basic approximation that 
we employ is that the fundamental form of the giant cusp solution is not seriously modified 
by the existence of an excess number of poles. Of course this approximation breaks down 
quantitatively already with one excess pole. Qualitatively however it holds well until the 
excess number of poles is of the order of the original number N{L) of the giant cusp solution. 
Another approximation is that the rest of the linear modes play no role in this case. At this 
point we hmit the discussion therefore to the situation 5N <S N{L) (regime II). 

To estimate the parameter A in the effective potential we consider the dynamics of one 
pole whose y position is far above ymax- According to Eq.(l.ll) the dynamics reads 

dya ^ 2v{N{L) + 5N) 1 

~d[ L ^^-^^^ 

Since the N{L) term cancels against the L"^ term (cf. Sec. II A), we remain with a repulsive 
term that in the effective potential translates to 

A = . (2.60) 

Next we estimate the value of the potential at the break-even point between attraction and 
repulsion. In the last subsection we saw that a foreign pole has to be inserted below ymax in 
order to be attracted towards the real axis. Now we need to push the new pole below the 
position of the existing pole whose index is iV(L) — 5N . This position is estimated as in Sec 
III C by employing the TFH distribution function (2.23). We find 

As before, this implies a threshold value of the amplitude of single pole solution A„iax sin 9 
which is obtained from equating A^ax = i^e~^*^. We thus find in the present case A^ax ^ 
u^{SN)'^/L'^. Using again a cubic representation for the effective potential we find a — 
2/{3iy'^SN) and 

Repeating the calculation of the escape rate over the potential barrier we find in the present 
case 

it:~^exp-^^(^^)V/^^ . (2.63) 
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For a given noise amplitude / there is always a value of L and v for which the escape rate 
is of 0(1) as long as 5N is not too large. When 5N increases the escape rate decreases, and 
eventually no additional poles can creep into the system. The typical number 5N for fixed 
values of the parameters is estimated from equating the argument in the exponent to unity 

5N^ [fL^'/v')^"' . (2.64) 

We can sec that 6N is strongly dependent on noise /, in contrast to regime III. Let us 
find the conditions of transition from regime II to III, where we see the saturation oi SN with 
respect to noise /. 

(i) Wc use the expression A^ax = 4z/e~^*^ for the amplitude of the pole solution that 
equals to cosii0t)^i^"-cos6> ' however, this is correct only for the large number ysN- When ysN < 1, 



a better approximation is Amax = 4^- From the equation (2.61) we find that the boundary 



value y^AT = 1 corresponds to bN ^ N{L)/2. 

(ii) We use the expression y^jq fa 2 In ^t^^j^ , but for a large value of 5N a better ap- 
proximation that can be found the same way is ysN ~ (-^) ~ ■k'^v{n{l)-6n) [12]. 
These expressions give us nearly Ithe same result for 5N f» N{L)/2. 

Prom (i) and (ii) we can make the following conclusions: 

(a) The transition from regime II to regime III generally occurs for 6N ^ N{L)/2; 

(b) Using the new expressions in (i) and (ii) for the amplitude Amax and ysN, we can 
determine the noise ^ in regime III by 

L ~ ^^^-^^ ~ ^ -J^^^^ ~ - {NiL)-6Ny ^'-^^^ 

This expression defines a very slight dependence of 6N on the noise / for SN > N{L)/2, 
which explains the noise saturation of 6N for regime III. 

(c) The form of the giant cusp solution is governed by the poles that are close to zero 
with respect to y. For the regime III, N{L)/2 poles that have positions y < ysN=N{L)/2 — 1 
remain at this position. This result explains why the giant cusp solution cannot be seriously 
modified for regime III. 

From eq. (2.64) by using the condition 

SN^N{L)/2 (2.66) 
the boundary noise between regimes II and III can be found as 

hr^u^ . (2.67) 
The basic equation describing pole dynamics follows 

dN 5N , , 

IF = ^ ■ P-^*^* 

where ^ is the number of poles that appear in the unit time in our system, SN is the 
excess number of poles, and T is the mean lifetime of a pole (between appearing and merging 
with the giant cusp). Using the result of numerical simulations for ^ and (2.66) we can find 
for T 
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'ig ~ TT" ~ const (2-71) 



T = ^ ~ uL'- . (2.69) 

dt 

Thus the hfetime is proportional to u and depends on the system size L very shghtly. 

Moreover, the hfetime of a pole is defined by the lifetime of the poles that are in a cusp. 
Prom the maximum point of the linear part of Eq.(2.1 ), we can find the mean character size 
(Fig.9( [28])) 

A™ ~ z/ (2.70) 
that defines the size of our cusps. The mean number of poles in a cusp 

X 
2z/ 

does not depend on L and u. The mean number of cusps is 

N,i,^ — ^-. (2.72) 
ribig V 

Let us assume that some cusp exists in the main minimum of the system. The lifetime of 
a pole in such a cusp is defined by three parts. 

(I) Time of the cusp formation. This time is proportional to the cusp size (with In- 
corrections) and the pole number in the cusp (from pole motion equations) 

Ti ~ Xmrihig ~ V (2.73) 

(II) Time that the cusp is in the minimum neighborhood. This time is defined by 

T2 - - (2.74) 

V 

where a is a neighborhood of minimum, such that the force from the giant cusp is smaller 
than the force from the fluctuations of the excess pole number 5N , and v is the velocity of a 
pole in this neighborhood. Fluctuations of excess pole number 5N are expressed as 

Nfi = Vm . (2.75) 
Prom this result and the pole motion equations we find that 

„^-N„^-^-^^-. (2.76) 

The velocity from the giant cusp is defined by 

v , , a a , , 

V ~ -N{L)- ~ - . (2.77) 
Ij Ij Ij 

So from equating these two equations we obtain 

a ~ VuL . (2.78) 

Thus for T2 we obtain 
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T2 ~ - ~ L . (2.79) 
(III) Time of attraction to the giant cusp. Prom the equations of motion for the poles we 

get 

Ta ~ Lln(-) - Lln\/L - L . (2.80) 
a 

The investigated domain of the system size was found to be 

Ti>T2,T3 (2.81) 

Therefore full lifetime is 

T = Ti + ^2 + ^3 ~ 1/ + sL , (2.82) 

where s is a constant and 

< s < 1 . (2.83) 

This result qualitatively and partly qualitatively explains dependence (2.69). From (2.69), 
(2.68), (2.66) we can see that in regime III ^ is saturated with the system size L. 

2.5.4 The acceleration of the flame front because of noise 

In this section we estimate the scaling exponents that characterize the velocity of the flame 
front as a function of the system size. To estimate the velocity of the flame front we need to 
create an equation for the mean oi < dh/dt > given an arbitrary number of poles in the 
system. This equation follows directly from (2.4) 

(dh\l_lr^' (2.84) 

\dt/ L^2t:Jo ^ ^ 

After substitution of (2.8) in (2.84) we get, using (2.11) and (2.12) 



Estimating the second and third terms in this equation are straightforward. Writing = 
N{L) + SN{L) and remembering that N{L) ^ L/u and 5N{L) ^ N{L)/2, we flnd that 
these terms contribute 0(1). The flrst term contributes only when the current of the poles 
is asymmetric. Noise introduces poles at a flnite value of i/min, whereas the rejected poles 
stream towards inflnity and disappear at the boundary of nonlinearity deflned by the position 
of the highest pole as 



2 In 



(2.86) 



Thus we have an asymmetry that contributes to the velocity of the front. To estimate 
the flrst term let us deflne 
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where E^''' ^ is the sum over the poles that are on the interval y : [1,1 + dl]. We can 
write 

^(E ^) = d{E + d{E ^ho^n , (2.88) 

where d{Y, ^)up is the flux of poles moving up and d{Y, ^)down is the flux of poles 
moving down. 

For these fluxes we can write 

din ^U, -rf(E ^)^n < ^dl . (2.89) 

So for the first term 



< E ^ = / ^%^<ii (2.90) 



/ 



k = l dt Jymin dl 



' Vmin dl 

^dN _ 

— yUmax ymin) 

dN 

— TrUmax 

dt 

Because of slight (In) dependence of ymax on L and v, ^ term determines order of 
nonlinearity for the first term in eq (2.85). This term equals zero for the symmetric current of 
poles and achieves the maximum for the maximal asymmetric current of poles. A comparison 
of V ~ L°-^V°°^ and ^ - L0.8fom confirms this calculation. 



2.6 Summary and Conclusions 

The main two messages of this chapter are: (i) There is an important interaction between the 
instability of developing fronts and random noise; (ii) This interaction and its implications 
can be understood qualitatively and sometimes quantitatively using the description in terms 
of complex poles. 

The pole description is natural in this context firstly because it provides an exact (and 
effective) representation of the steady state without noise. Once one succeeds to describe 

also the perturbafAons about this steady state in terms of poles, one achieves a particularly 
transparent language for the study of the interplay between noise and instability. This 
language also allows us to describe in qualitative and semi-quantitative terms the inverse 
cascade process of increasing typical lengths when the system relaxes to the steady state 
from small, random initial conditions. 
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The main conceptual steps in this chapter are as follows: firstly one realizes that the 
steady state solution, which is characterized by N{L) poles aligned along the imaginary axis 
is marginally stable against noise in a periodic array of L values. For all values of L the 
steady state is nonlinearly unstable against noise. The main and foremost effect of noise 
of a given amphtude / is to introduce an excess number of poles SN{L, /) into the system. 
The existence of this excess number of poles is responsible for the additional wrinkling of the 
flame front on top of the giant cusp, and for the observed acceleration of the flame front. By 
considering the noisy appearance of new poles we rationalize the observed scaling laws as a 
function of the noise amplitude and the system size. 

Theoretically we therefore concentrate on estimating 5N(L, /). We note that some of our 
consideration are only qualitative. For example, we estimated SN{L, f) by assuming that 
the giant cusp solution is not seriously perturbed. On the other hand we flnd a flux of poles 
going to infinity due to the introduction of poles at finite values of y by the noise. The 
existence of poles spread between i/max and infinity is a significant perturbation of the giant 
cusp solution. Thus also the comparison between the various scaling exponents measured 
and predicted must be done with caution; we cannot guarantee that those cases in which our 
prediction hit close to the measurement mean that the theory is quantitative. However we 
beheve that our consideration extract the essential ingredients of a correct theory. 

The "phase diagram" as a function of L and / in this system consists of three regimes. 
In the first one, discussed in Section 2.5.3 , the noise is too small to have any effect on the 
giant cusp solution. In the second the noise introduces excess poles that serve to decorate the 
giant cusp with side cusps. In this regime we find scaling laws for the velocity as a function 
of L and / and we are reasonably successful in understanding the scaling exponents. In the 
third regime the noise is large enough to create small scale structures that are not neatly 
understood in terms of individual poles. It appears from our numerics that in this regime 
the roughening of the flame front gains a contribution from the the small scale structure in 
a way that is reminiscent of stable, noise driven growth models like the Kardar-Parisi- Zhang 
model. 

One of our main motivations in this research was to understand the phenomena observed 
in radial geometry with expanding flame fronts. . We note that many of the insights offered 
above translate immediately to that problem. Indeed, in radial geometry the flame front 
accelerates and cusps multiply and form a hierarchic structure as time progresses. Since the 
radius (and the typical scale) increase in this system all the time, new poles will be added 
to the system even by a vanishingly small noise. The marginal stability found above holds 
also in this case, and the system will allow the introduction of excess poles as a result of 
noise. The results discussed in Ref. [19] can be combined with the present insights to provide 
a theory of radial growth ( chapter 4). 

Finally, the success of this approach in the case of flame propagation raises hope that 
Laplacian growth patterns may be dealt with using similar ideas. A problem of immediate 
interest is Laplacian growth in channels, in which a flnger steady-state solution is known 
to exist. It is documented that the stability of such a flnger solution to noise decreases 
rapidly with increasing the channel width. In addition, it is understood that noise brings 
about additional geometric features on top of the flnger. There are enough similarities here 
to indicate that a careful analysis of the analytic theory may shed as much light on that 
problem as on the present one. 
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Chapter 3 



Using of Pole Dynamics for Stability 
Analysis of Flame Fronts: Dynamical 
Systems Approach in the Complex 
Plane 



3.1 Introduction 

In this chapter we discuss the stabihty of steady flame fronts in channel geometry. We write 
shortly about this topic in chapter 2 (Sec. 2.3) and we want to consider it in detail in this 
chapter. Traditionally [1-3] one studies stability by considering the linear operator which is 
obtained by linearizing the equations of motion around the steady solution. The eigenfunc- 
tions obtained are delocalized and in certain cases are not easy to interpret. In the case of 
flame fronts the steady state solution is space dependent and therefore the eigenfunctions are 
very different from simple Fourier modes. We show in this chapter that a good understand- 
ing of the nature of the eigenspectrum and eigenmodes can be obtained by doing almost the 
opposite of traditional stability analysis, i.e., studying the localized dynamics of singularities 
in the complex plane. By reducing the stability analysis to a study of a finite dimensional 
dynamical system one can gain considerable intuitive understanding of the nature of the 
stability problem. 

The analysis is based on the understanding that for a given channel width L the steady 
state solution for the flame front is given in terms of N{L) poles that are organized on a line 
parallel to the imaginary axis [12]. Stability of this solution can then be considered in two 
steps. In the first step we examine the response of this set of N{L) poles to perturbations 
in their positions. This procedure yields an important part of the stability spectrum. In the 
second step we examine general perturbations, which can also be described by the addition 
of extra poles to the system of N(L) poles. The response to these perturbations gives us 
the rest of the stability spectrum; the combinations of these two steps rationalizes all the 
qualitative features found by traditional stability analysis. 

In Sec. 2 we present the results of traditional linear stability analysis, and show the eigen- 
values and eigenfunctions that we want to interpret by using the pole decomposition. Sec. 
3 presents the analysis in terms of complex singularities, in two steps as discussed above. A 
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summary and discussion is presented in Sec.4. 



3.2 Linear Stability Analysis in Channel Geometry 

The standard technique to study the hnear stabihty of the steady solution is to perturb it 
by a small perturbation 0(6', t): u{9, t) = UgiO) + (j){6, t) . Linearizing the dynamics for small 
results in the following equation of motion 



d<i){e,t) _ 1 

dt ~ 1? 



de[us{e)(t>{e,t)\ 



+ udl4>{e,t)\ + ^i{4>{e,t)) . (3.1) 

were the linear operator contains Us{0) as a coefficient. Accordingly simple Fourier modes do 
not diagonalize it. Nevertheless, we proceed to decompose 0(x) in Fourier modes according 
to , 

0(^,t) = f; 4^k{t)e''' (3.2) 

k=—oo 

oo N 

Us{e) = -2ui ^si^n(A;)e-l'=l^^e''=^ (3.3) 

fc=— oo j=l 

The last equation follows from (2.13) by expanding in a series of sin k6. In these sums the 
discrete k values run over all the integers. Substituting in Eq.(3.1) we get: 

^^-Y^knUt), (3.4) 

where akn are entires of an infinite matrix: 

akk = - ' (3-5) 

k ^ 
akn = -^sign{k-n){2i^J2^^^''~''^'") ^^n. (3.6) 

To solve for the eigenvalues of this matrix we need to truncate it at some cutoff A;- vector k*. 
The scale k* can be chosen on the basis of Eq.(3.5) from which we see that the largest value 
of k for which a^^ > is a scale that we denote as kmaxi which is the integer part oiLjv. We 
must choose k* > kmax and test the choice by the convergence of the eigenvalues. The chosen 
value of k* in our numerics was Ak^ax ■ One should notice that this cutoff limits the number of 
eigenvalues, which should be infinite. However the lower eigenvalues will be well represented. 
The results for the low order eigenvalues of the matrix akn that were obtained from the 
converged numerical calculation are presented in Fig.3.1 The eigenvalues are multiplied by 
L'^/u and are plotted as a function of L. We order the eigenvalues in decreasing order and 
denote them as Aq > Ai > A2 In addition to the eigenvalues, the truncated matrix also 
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Figure 3.1: A plot of the first five eigenvalues obtained by diagonalizing the matrix obtained 
by traditional stability analysis, against the system size. The eigenvalues are normalized by 
1? jv. The largest eigenvalue is zero, which is a Goldstone mode. All the other eigenvalues are 
negative except for the second and third that touch zero periodically. The second and fourth 
eigenvalues are represented by a solid line and the third and fifth eigenvalues are represented 
by a dot-dashed fine. 

yields eigenvectors that we denote as AS^^ . Each such vector has k* entries, and we can 
compute the eigenfunctions of the linear operator (3.1), using (3.2), as 

f\B)^Y.^^'Af . (3.7) 

-k* 

Eq.(3.1) docs not mix even with odd solutions in as can be checked by inspection. Conse- 
quently the available solutions have even or odd parity, expandable in cither cos or sin func- 
tions. The first two nontrivial eigenfunctions f^^^{9) and f^'^\6) are shown in Figs. 3. 2, 3. 3. 
a It is evident that the function in Fig. 3. 2 is odd around zero whereas in Fig.3.3 it is even. 
Similarly we can numerically generate any other eigenfunction of the linear operator, but we 
understand neither the physical significance of these eigenfunction nor the L dependence of 
their associated eigenvalues shown in Fig.3.1 In the next section we will demonstrate how the 
dynamical system approach in terms of singularities in the complex plane provides us with 
considerable intuition about these issues. 



3.3 Linear Stability in terms of complex singularities 

Since the partial differential equation is continuous there is an infinite number of modes. To 
understand this in terms of pole dynamics we consider the problem in two steps: First, we 
consider the 2N{L) modes associated with the dynamics of the N{L) poles of the giant cusp. 
In the second step we explain that all the additional modes result from the introduction of 
additional poles, including the reaction of the N{L) poles of the giant cusp to the new poles. 
After these two steps we will be able to identify all the linear modes that were found by 
diagonalizing the stability matrix in the previous section. 
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Figure 3.2: The first odd eigenfunction obtained from traditional stability analysis. 
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Figure 3.3: The first even eigenfunction obtained from traditional stability analysis. 
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3.3.1 The modes associated with the giant cusp 



In the steady solution all the poles occupy stable equibirilium positions. The forces operating 
on any given pole cancel exactly, and we can write matrix equations for small perturbations 
in the pole positions 5yi and 5xi. 

Following [12] we rewrite the equations of motion (2.12) using the Lyapunov function U: 

dU 

Lm = ^ (3.8) 
dyi 

where i — 1, N and 



^[^lnsinh|/i + 2^(lnsinh^^^ 
^ i i<k 2 



+ lnsinh^^)]-5:i/, (3.9) 

i 

The linearized equations of motion for Syi are: 

The matrix d'^U/dyidyk is real and symmetric of rank N. We thus expect to find N real 
eigenvalues and N orthogonal eigenvectors. 

For the deviations 6xi in the x positions we find the following linearized equations of 
motion 



N ^ 



L ^^f^^. cosh{yj -yk)-l) 
1 



In shorthand: 



cosh(|/j + i/fc) - 1 

u ^ 1 1 

L kJ^j^j cosh{yj - yu) - 1) cosh(yj + y^) - 1 

dox ■ 

L-^ = VikSxk . (3.12) 

The matrix V is also real and symmetric. Thus V and S^U / dyidyk together supply 2N{L) 
real eigenvalues and 2N{L) orthogonal eigenvectors. The explicit form of the matrices V and 
d'^U / dyidyk is as follows: For i k: 

_^.^Z2 IJ^ 



dyidyk L^smh\^) sinh2(^) 

= 7( U \ 1^ + 1 ^ J (3-14) 

L cosh(7/i - yk) - 1) cosh(yj + y^) - 1 



47 




Figure 3.4: The eigenvalues associated with perturbing the positions of the poles that consist 
the giant cusp. The largest eigenvalue is zero. The second, third, fourth and fifth eigenvalues 
are represented by a solid line, dot-dashed line, dotted line and dashed line respectively. 



and ior i — k one gets: 




+ 



dyj L^;^. V2sinh'(Mi) 2s\nh^{y^)) 

^ sinh^(yj)^ 

^ u 1 1 

= 5^"Z^cosh(y,-y,)-l) + cosh(y, + y,)-l^^ ^^'^^^ 

Using the known steady state solutions yi at any given L we can diagonalize the N{L) x N{L) 
matrices numerically. In Fig.3.4 we present the eigenvalues of the lowest order modes obtained 
from this procedure. The least negative eigenvalues touch zero periodically. This eigenvalue 
can be fully identified with the motion of the highest pole yN{L) in the giant cusp. At isolated 
values of L the position of this pole tends to infinity, and then the row and the column in our 
matrices that contain yN(L) vanish identically, leading to a zero eigenvalue. The rest of the 
upper eigenvalues match perfectly with half of the observed eigenvalues in Fig. 3.1. In other 
words, the eigenvalues observed here agree perfectly with the ones plotted in this Fig. 3.1 until 
the discontinuous increase from their minimal points. The "second half" of the oscillation in 
the eigenvalues as a function of L is not contained in this spectrum of the N(L) poles of the 
giant cusp. To understand the rest of the spectrum we need to consider perturbation of the 
giant cusp by additional poles. The eigenfunctions can be found using the knowledge of the 
eigenvectors of these matrices. Let us denote the eigenvectors of d^U/dyidyk and V as a^^^ 
and b^^^ respectively. The perturbed solution is explicitly given as (taken for Xs = 0): 

Us{0) + Su 2z/g ^^^^^^^ ^ ^^^^ _ ^^^^^ _ (3.17) 

where 5u is 

N oo 

Su^ - 4z/ ^ ^ 5yike~^'"' sin kO 

i=l k=l 
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Figure 3.5: The first odd eigenfunction associated with perturbing the positions of the poles 
in the giant cusp. 

N oo 

- 4z/^^5xiA;e-^J''cosA;^ (3.18) 

i=l k=l 

So knowing the eigenvectors a'-^^ and fe*^^-* we can estimate the eigenvectors f'^^\9) of (3.7): 

N oo 

/i£(^^) = -4z/^5:a?)A;e-^^'sinfc^, j = l,...,iV (3.19) 

i=l k=l 

or 

N oo 

f^^{e) = -AuY.Y.^^pke-^y^coske , j = i,...,iv (3.20) 

i=l fc=l 

where we display separately the sin expansion and the cos expansion. For the case j = 1, 
the eigenvalue is zero, and a uniform translation of the poles in any amount 5xi results in a 
Goldstone mode. This is characterized by an eigenvector h^P = 1 for all i. The eigenvectors 
(Fig.3.5,3.6)computed this way are identical to numerical precision with those shown in 
Figs. 3. 2, 3. 3, and observe the agreement. 

3.3.2 Modes related to additional poles 

In this subsection we identify the rest of the modes that were not found in the previous 
subsection. To this aim we study the response of the TFH solution to the introduction 
of additional poles. We choose to add M new poles all positioned at the same imaginary 
coordinate i/p -C Umax, distributed at equidistant real positions {xj = Xq + (27r/M)j}^^. For 
Xo = we use (2.8) and the Fourier expansion to obtain a perturbation of the form 

6u{e, t) ~ AvMe~^"^y^^*^ sin MO (3.21) 

For xq = —tt/2M we get 

5u{e, t) ~ AvMe-^^y^^*^ cos MO (3.22) 
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Figure 3.6: The first even eigenfunction associated with perturbing the positions of the poles 
in the giant cusp. 



in both cases the equations for the dynamics of Up follow from Eqs.(2.11)-(2.12): 

^ ^ 2-^a(M) , (3.23) 

where a{M) is given as: 

a(M) = [i(^ + l)]-i(^-M) (3.24) 

Since (3.23) is linear, we can solve it and substitute in Eqs.(3.21)-(3.22). Seeking a form 
6u{9,t) ~ exp{—\{M)t) we find that the eigenvalue A(M) is 

A(M) = 2M-^a(M) (3.25) 

These eigenvalues are plotted in Fig. 3. 7 At this point we consider the dynamics of the poles 
in the giant cusp under the influence of the additional M poles. From Eqs.(3.10), (3.12), 
(2.11), (2.12) we obtain, after some obvious algebra, 

L^y^ = E if^^Vi - 4^Me-*^^''W sinh(M|/,) (3.26) 

or 

LSxi = VijSxj - 4yMe-^^''W cosh{Myi) (3.27) 

It is convenient now to transform from the basis 6yi to the natural basis Wi which is obtained 
using the linear transformation w = A~^6y. Here the matrix A has columns which are 
the eigenvectors of d^U/dyidyj which were computed before. Since the matrix was real 
symmetric, the matrix A is orthogonal, and A^'^ = A'^. Define C = 4-^Me~*^^''''°^ and write 



Wi 



-Kw, - Ce-"(*^)*ei , (3.28) 
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Figure 3.7: Spectrum of eigenvalues associated with the reaction of the poles in the giant 
cusp to the addition of new poles. 



where — are the eigenvalues associated with the columns of A, and 



(3.29) 



We are looking now for a solution that decays exponentially at the rate A(M): 

Wi{t) = «;i(0)e-^^^)* (3.30) 
Substituting the desired solution in (3.28) we find a condition on the initial value of Wi'. 

C 



Xi - A(M) 



(3.31) 



Transforming back to Syi we get 



C 



k k Xk-A{M) I 



Z Aik sinh Myi 



-C^sinhMy^Z 



AikAik 



Afe - A^ 



(3.32) 



We can get the eigenfunctions of the hnear operator, as before, using Eqs.(3.18), (3.21), 
(3.22), (3.32). We get 



NiL) 



A^rnA 



im-^lm 



X ke-^^' sin kO + L^C sin MO 



i=l k=l 
2/ 



A^ - A(M) 



) 



(3.33) 



An identical calculation to the one started with Eq. (3.28) can be followed for the deviations 
5xi. The final result reads 



N{L) oo 



A ■ A 



1=1 k=l I m Xm — A(M) 

X ke-'^y^ cos ke + L'^C cos MO , 



(3.34) 
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where A is the matrix whose columns are the eigenvectors of V, and — Aj its eigenvalues. 

We are now in position to explain the entire linear spectrum using the knowledge that we 
have gained. The spectrum consists of two separate types of contributions. The first type 
has 2N modes that belong to the dynamics of the unperturbed N{L) poles in the giant cusp. 
The second part, which is most of the spectrum, is built from modes of the second type since 
M can go to infinity. This structure is seen in the Fig. 3. 4 and Fig.3.7. 

We can argue that the set of eigenfunctions obtained above is complete and exhaustive. 
To do this we show that any arbitrary periodic function of 6 can be expanded in terms of 
these eigenfunctions. Start with the standard Fourier series in terms of sin and cos functions. 
At this point solve for sin k9 and cos k9 from Eqs. (3.33-3.34). Substitute the results in the 
Fourier sums. We now have an expansion in terms of the eigenmodes /(*^) and in terms of 
the triple sums. The triple sums however can be expanded, using Eqs. (3. 19-3.20), in terms of 
the eigenfunctions f^^\ We can thus decompose any function in terms of the eigenfunctions 
/(M) and /W. 

3.4 Conclusions 

We discussed the stability of flame fronts in channel geometry using the representation of 
the solutions in terms of singularities in the complex plane. In this language the stationary 
solution, which is a giant cusp in configuration space, is represented by N{L) poles which 
are organized on a line parallel to the imaginary axis. We showed that the stability problem 
can be understood in terms of two types of perturbations. The first type is a perturbation 
in the positions of the poles that make up the giant cusp. The longitudinal motions of the 
poles give rise to odd modes, whereas the transverse motions to even modes. The eigenvalues 
associated with these modes are eigenvalues of a finite, real and symmetric matrices, cf. 
Eqs. (3. 13), (3.14), (3.15), (3.16). The second type of perturbations is obtained by adding 
poles to the set of N{L) poles representing the giant cusp. The reaction of the latter poles is 
again separated into odd and even functions as can be seen from Eqs. (3. 21), (3.22). Together 
the two types of perturbations rationalize and explain all the features of the eigenvalues and 
eigenfunctions obtained from the standard linear stability analysis. 
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Chapter 4 

Dynamics and Wrinkling of Radially 
Propagating Fronts Inferred from 
Scaling Laws in Channel Geometries 

4.1 Introduction 



The main idea of this chapter is that in order to derive scahng laws for unstable front 
propagation in radial geometry, it is useful to study noisy propagation in channel geometries, 
in which the noiseless dynamics results usually in simple shapes of the advancing fronts [1]. 

The understanding of radial geometries requires control of the effect of noise on the 
unstable dynamics of propagation. It is particularly difficult to achieve such a control in 
radial geometries due to the vagueness of the distinction between external noise and noisy 
initial conditions. Channel geometries are simpler when they exhibit a stable solution for 
growth in the noiseless limit. One can then study the effects of external noise in such 
geometries without any ambiguity. If one finds rules to translate the resulting understanding 
of the effects of noise in channel growth to radial geometries, one can derive the scaling 
laws in the later situation in a satisfactory manner. We will exemplify the details of such a 
translation in the context of premixed flames that exist as self sustaining fronts of exothermic 
chemical reactions in gaseous combustion. But our contention is that similar ideas should 
be fruitful also in other contexts of unstable front propagation. Needless to say, there are 
aspects of the front dynamics and statistics in the radial geometry that cannot be explained 
from observations of fronts in a channel geometry; examples of such aspects are discussed at 
the end of this chapter. 

Mathematically our example is described [11] by an equation of motion for the angle- 
dependent modulus of the radius vector of the flame front, R{9, t): 



Here 0< ^ < 27r is an angle and the constants Ub, Dm and 7 are the front velocity for an ideal 
cylindrical front, the Markstein diffusivity and the thermal expansion coefficient respectively. 



dR 
'dt 



+ 




(4.1) 
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Ro{t) is the mean radius of the propagating flame: 

Ro{t)^^ r R{e,t)de . (4.2) 

ZTT JO 

The functional I{R) is best represented in terms of its Fourier decomposition. Its Fourier 
component is where is the Fourier component of R. Simulations of this equation, as 
well as experiments in the parameter regime for which this equation is purportedly relevant, 
indicate that for large times Rq grows as a power in time 

Ro{t) = [const + tf , (4.3) 

with /? > 1, and that the width of the interface W grows with Rq as 

W{t) ~ R^{tY , (4.4) 

with X < 1. 



4.2 The Geometry of Developing Flame Fronts: Anal- 
ysis with Pole Decomposition 

The study of growing fronts in nonlinear physics [1] offers fascinating examples of spontaneous 
generation of fractal geometry [2,3]. Advancing fronts rarely remain flat; usually they form 
either fractal objects with contorted and ramified appearance, like Laplacian growth patterns 
and diffusion limited agregates (DLA) [31], or they remain graphs, but they "roughen" in 
the sense of producing self-affine fractals whose "width" diverges with the linear scale of 
the system with some characteristic exponent. The study of interface growth where the 
roughening is caused by the noisy environment, with either annealed or quenched noise, was 
a subject of active research in recent years [32,33]. These studies met considerable success 
and there is significant analytic understanding of the nature of the universality classes that 
can be expected. The study of interface roughening in system in which the flat surface is 
inherently unstable is less developed. One interesting example that attracted attention is 
the Kuramoto-Sivashinsky equation [9, 34] which is known to roughen in 1+1 dimensions 
but is claimed not to roughen in higher dimensions [42]. Another outstanding example is 
Laplacian growth patterns [35] . This chapter is motivated by a new example of the dynamics 
of outward propagating flames whose front wrinkles and fractalizes [11]. We will see that this 
problem has many features that closely resemble Laplacian growth, including the existence of 
a single flnger in channel growth versus tip splitting in cylindrical outward growth, extreme 
sensitivity to noise, etc. In the case of flame fronts the equation of motion is amenable to 
analytic solutions and as a result we can understand some of these issues. 

The physical problem that motivates this analysis is that of pre-mixed flames which exist 
as self-sustaining fronts of exothermic chemical reactions in gaseous combustion. It had been 
known for some time that such flames are intrinsically unstable [43]. It was reported that 
such flames develop characteristic structures which includes cusps, and that under usual 
experimental conditions the flame front accelerates as time goes on [10]. In recent work 
Filyand et al. [11] proposed an equation of motion that is motivated by the physics and 
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seems to capture a number of the essential features of the observations. The equation is 
written in cyhndrical geometry and is for R{6, t) which is the modulus of the radius vector 
on the flame front: 

dR ^ Ut (OrV Dm d^R , . 

dt 27?o'W [dej Ro\t) 86^ ^ ' ' 

Here 0< ^ < 27r is an angle and the constants Ub, Dm and 7 are the front velocity for an ideal 
cylindrical front, the Markstein diffusivity and the thermal expansion coefficient respectively. 
Ro{t) is the mean radius of the propagating flame: 

Ro{t) = ^ r R{e,t)de . (4.6) 

ZTT Jo 

The functional I{R) is best represented in terms of its Fourier decomposition. Its Fourier 
component is \k\Rk where Rk is the Fourier component of R. 

Numerical simulations of the type reported in ref. [11] are presented in Fig.4.1. The two 
most prominent features of these simulations are the wrinkled multi-cusp appearance of the 
fronts and its acceleration as time progresses. One observes the phenomenon of tip splitting 
in which new cusps arc added to the growing fronts between existing cusps. Both experiments 
and simulations indicate that for large times Rq grows as a power in time 

Ro{t) = {const + tf , (4.7) 

with /3 > 1, (of the order of 1.5) and that the width of the interface W grows with Rq as 

W{t) ~ i?o(i)^ , (4.8) 

with X < 1 (of the order of 2/3). The understanding of these two features and the derivation 
of the scaling relation between (3 and x the main aims of this chapter. 

Equation (4.6) can be written as a one-parameter equation by rescaling R and t according 
to r = RUb/DM, T = HJ^/Dm- Computing the derivative of Eq.(4.6) with respect to 9 and 
substituting the dimensionless variables one obtains: 

or Tq do Tq oO^ 2ro 

where u = To complete this equation we need a second one for ro(t), which is obtained 
by averaging (4.6) over the angles and rescaling as above. The result is 



dr 2rQ 27r 7o 
These two equations are the basis for further analysis 
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Following [12,14-16,38,39] we expand now the solutions u{9,t) in poles whose position 
Zj{T) = Xj{T) + iyjij) in the complex plane is time dependent: 



N 




u{e,T) 



Ecot 



+ c.c. 



(4.11) 



2sm[e-Xj{T)] 



p{ cosh[|/j(r)] - cos[^ - Xj{T)] ' 



N 



r{e,T) = 2^1n [cosh(|/^(r)) - cos(^ - Xjir))] + C(t) . 



(4.12) 



In (4.12) C(t) is a function of time. The function (4.12) is a superposition of quasi-cusps 
(i.e. cusps that arc rounded at the tip). The real part of the pole position (i.e. xj) describes 
the angle coordinate of the maximum of the quasi-cusp, and the imaginary part of the pole 
position (i.e yj) is related the height of the quasi-cusp. As yj decreases (increases) the height 
of the cusp increases (decreases). The physical motivation for this representation of the 
solutions should be evident from Fig.4.1. 

The main advantage of this representation is that the propagation and wrinkling of the 
front can be described now via the dynamics of the poles and of ro(t). Substituting (4.11) in 
(4.9) we derive the following ordinary differential equations for the positions of the poles: 



After substitution of (4.11) in (4.10) we get, using (4.13) the ordinary differential equation 
for ro, 



In our problem the outward growth introduces important modifications to the channel 

results. The number of poles in a stable configuration is proportional here to the radius 
Tq instead of L, but the former grows in time. The system becomes therefore unstable to 
the addition of new poles. If there is noise in the system that can generate new poles, they 
will not be pushed toward infinite y. It is important to stress that any infinitesimal noise 
(either numerical or experimental) is sufficient to generate new poles. These new poles do 
not necessarily merge their x-positions with existing cusps. Even though there is attraction 
along the real axis as in the channel case, there is a stretching of the distance between the 
poles due to the radial growth. This may counterbalance the attraction. Our first new idea 
is that these two opposing tendencies define a typical scale denoted as C if we have a cusp 
that is made from the a;-merging of A^^ poles on the line x = Xc and we want to know whether 
a x-nearby pole with real coordinate Xi will merge with this large cusp, the answer depends 
on the distance D = ro\xc — Xi\. There is a length C{Nc, Tq) such that if D > C{Nc, ro) then 
the single cusp will never merge with the larger cusp. In the opposite limit the single cusp 
will move towards the large cusp until their x- position merges and the large cusp will have 
Nc + 1 poles. 




27V 



(4.13) 




(4.14) 
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This finding stems directly from tlie equations of motion of the x-merged poles and 
the single pole at Xi. First note that from Eq.4.7 (which is not explained yet) it follows that 
asymptotically ro(T) = (a + r)'^ where ro(0) = a^. Next start from 4.13 and write equations 
for the angular distance x — Xi — Xc- It follows that for any configuration yj along the 
imaginary axis 

dx 2Ar^sinx[l - cosx]-^ _ 2iVcCot(f) 

d^- (a + rfP ~ ~ (a + rf/^ ' ^ ' 



For small x we get 



The solution of this equation is 



To find C we set a;(oo) > from which we find that the angular distance will remain finite 
as long as 

^(0)^ > ■ (4-18) 

Since vq ~ we find the threshold angle x* 

X* ~ ViVcro , (4.19) 

above which there is no merging between the giant cusp and the isolated pole. To find the 
actual distance jC{Nc, ro) we multiply the angular distance by ro and find 



£{N„ ro) = rox* - ^jNeV^' . (4.20) 

To understand the geometric meaning of this result we recall the features of the TFH cusp 
solution. Having a typical length L the number of poles in the cusp is linear in L. Similarly, 
if we have in this problem two cusps a distance 2£ apart, the number in each of them 
will be of the order of C From (4.20) it follows that 

Cr^4 . (4.21) 

For /3 > 1 the circumference grows faster than C, and therefore at some points in time poles 
that appear between two large cusps would not be attracted toward either, and new cusps 
will appear. We will show later that the most unstable positions to the appearance of new 
cusps are precisely the midpoints between existing cusps. This is the mechanism for the 
addition of cusps in analogy with tip splitting in Laplacian growth. 

We can now estimate the width of the flame front as the height of the largest cusps. Since 
this height is proportional to C (cf. property (v) of the TFH solution), Eq.(4.21) and Eq.(4.8) 
lead to the scaling relation 

X = l/(3. (4.22) 

This scaling law is expected to hold all the way to (3 = 1 for which the flame front does not 
accelerate and the size of the cusps becomes proportional to tq. 
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In channels there is a natural lengthscale, the width L of the channel. The translation of 
channel results to radial geometry will be based on the identification in the latter context of 
the time dependent scale C{t) that plays the role of L in the former. To do this we need first 
to review briefly the main pertinent results for noisy channel growth. In channel geometry 
the equation of motion is written in terms of the position h{x^ t) of the flame front above the 
X-axis. After appropriate rescalings [12] it reads: 



dh{x,t) _ 1 



dh{x, t) 
dx 



+u^^^ + I{h{x,t)} + l. (4.23) 



It is convenient to rescale the domain size further to < ^ < 27r, and to change variables to 
u{9, t) = dh{9, t)/d6. In terms of this function we find 

du{o,t) u(e,t) du{e,t) V dMo,t) i 

where L — L/2t:. In noiseless conditions this equation admits exact solutions that are 
represented in terms of N poles whose position Zj{t) = Xj{t) + iyj{t) in the complex plane is 
time dependent: 

^ ^9-Zj{ty 



u{9, t) = cot 



=1 



+ c.c. , (4.25) 



The steady state for channel propagation is unique and linearly stable; it consists of N{L) 
poles which are aligned on one line parallel to the imaginary axis. The geometric appearance 
of the flame front is a giant cusp, analogous to the single finger in the case of Laplacian growth 
in a channel. The height of the cusp is proportional to L, and the propagation velocity is a 
constant of the motion. The number of poles in the giant cusp is linear in L Eq. (1.10, 

The introduction of additive random noise to the dynamics changes the picture quali- 
tatively. It is convenient to add noise to the equation of motion in Fourier representation 
by adding a white noise rjk for every k mode. The noise correlation function satisfies the 
relation < r)k{t)r)k'{t') >— Sk^k'^{t — t')f/L. The noise in our simulations is taken from a 



flat distribution in the interval [—y2f/L,y2f/L]; this guarantees that when the system size 
changes, the typical noise per unit length of the fiame front remains constant. It was shown 
in chapter 2 and [17] that for moderate but fixed noise levels the average velocity v of the 
front increases with L as a power law. In our present simulations we found 

vr^ L" , // 0.35 ± 0.03 . (4.26) 

For a fixed system size L the velocity has also a power law dependence on the level of the 
noise, but with a much smaller exponent: v ~ , 0.02. These results were understood 

theoretically by analysing the noisy creation of new poles that interact with the poles defining 
the giant cusp (in chapter 2 and [17]). 

Next we shed light on phenomenon of tip sphtting that here is seen as the addition of new 
cusps roughly in between existing ones. We mentioned the instability toward the addition of 
new poles. We argue now that the tip between the cusps is most sensitive to pole creation. 
This can be shown in both channel and radial geometry. For example consider a TFH-giant 
cusp solution in which all the poles are aligned (without loss of generality) on the x = line. 
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Add a new pole in the complex position {xa, Ua) to the existing N{L) poles, and study its 
fate. It can be shown that in the limit Ua ^ oo (which is the limit of a vanishing perturbation 
of the solution) the equation of motion is 

^ = ^(2iV(L) + 1) - 1 I/. ^ oo . (4.27) 

Since N{L) satisfies (1.10) this equation can be rewritten as 

dya 47rz/ 



dt L2 



1-a) Va^oo, (4.28) 



where a = (L/(27rt/) + l)/2 — N{L). Obviously a < 1 and it is precisely 1 only when L is 
L — {2n + 1)2^1/. Next it can be shown that for much larger than yN{L) but not infinite 
the following is true: 

^ > hm ^ x„ = (4.29) 

dya ^ y dya , . 

— < hm — = TT (4.30) 

dt Va^OO dt 

We learn from these results that there exist values of L for which a pole that is added at 
infinity will have marginal attraction {dya/dt = 0). Similar understanding can be obtained 
from a standard stability analysis without using pole decomposition. Perturbing a TFH- 
cusp solution we find linear equations whose eigenvalues Aj can be obtained by standard 
numerical techniques: (i) all Re{Xi) are non-positive, (ii) at the isolated values of L for 
which L = (2n + l)27rz/ i?e(Ai) and Re{X2) become zero (note that due to the logarithmic 
scale the zero is not evident) (iii) There exists a general tendency of all Re{\i) to approach 
zero in absolute magnitude such as from below as L increases. This indicates a growing 
sensitivity to noise when the system size increases, (iv) There exists a Goldsone mode Aq = 
due to translational invariance. 

The upshot of this discussion is that finite perturbations (i.e. poles at finite ya) will grow 
if the X position of the pole is sufficiently near the tip. The position x = tt (the tip of the 
finger) is the most unstable one. In the channel geometry this means that noise results in 
the appearance of new cusps at the tip of the fingers, but due to the attraction to the giant 
cusp they move toward x = and disappear in the giant cusp. In fact, one sees in numerical 
simulations a train of small cusps that move toward the giant cusp. Analysis shows that 
at the same time the furthest pole at yN{L) is pushed towards infinity. Also in cylindrical 
geometry the most sensitive position to the appearance of new cusps is right between two 
existing cusps independently if the system is marginal (the total number of poles fits the 
radius) or unstable (total number of poles is too small at a given radius). Whether or not 
the addition of a new pole results in tip splitting depends on their x position. When the 
distance from existing cusps is larger than C the new poles that arc generated by noise will 
remain near the tip between the two cusps and will cause tip splitting. 

The picture used remains valid as long as the poles that are introduced by the noisy 
perturbation do not destroy the identity of the giant cusp. Indeed, the numerical simulations 
show that in the presence of moderate noise the additional poles appear as smaller cusps that 
are constantly running towards the giant cusp. Our point here, is not to predict the numerical 
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Figure 4.1: Simulations of the outward propagating flame front. Note that there is a wide 
distribution of cusp sizes. 

values of the scaling exponents in the channel (this was done in ref. [17] and chapter 2, but to 
use them to predict the scaling exponents characterizing the acceleration and the geometry 
of the flame front in radial geometry. 

Superficially it seems that in radial geometry the growth pattern is qualitatively different. 
In fact, close observation of the growth patterns (see Fig. 4.1) shows that most of the time there 
exist some big cusps that attract other smaller cusps, but that every now and then "new" 
big cusps form and begin to act as local absorbers of small cusps that appear randomly. The 
understanding of this phenomenon gives the clue how to translate results from channels to 
radial growth. 

Equations (4.9), (4.10) again admit exact solutions in terms of poles, of the form of 
Eq.(4.13). It is easy write down the equations of motion of the poles and check that the 
poles are attractive along the real direction (which means physically that they are attracted 
along the angular coordinate) but they are repulsive along the imaginary direction, which is 
associated with the radial coordinate. If it were not for the stretching that is caused by the 
increase of the radius (and with it the perimeter), all the poles would have coalesced into one 
giant cusp. Thus we have a competition between pole attraction and stretching. Since the 
attraction decreases with the distance between the poles in the angular directions, there is 
always an initial critical length scale above which poles cannot coalesce their real coordinates 
when time progresses. 

Suppose now that noise adds new poles to the system. The poles do not necessarily merge 
their real positions with existing cusps. If we have a large cusp made from the merging 
of the real coordinates Xc of poles, wc want to know whether a nearby pole with real 
coordinate Xi will merge with this large cusp. The answer will depend of course on the 
distance D = ro\xc — xi\. A direct calculation [19], using the equation of motions for the 
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Figure 4.2: Panel a: a logarithmic plot of the velocity versus time for a radially evolving 
system. The parameters of the simulation are: / = 10~^, 7 = 0.8, u = 1. Panel b:Logarithmic 
plot of the width of the flame front as a function of the mean radius. 

poles shows that there exists a critical length C{ro) such that if D > C{ro) the single pole 
never merges with the giant cusp. The result of the calculation is that 

C ~ rl^^ . (4.31) 

Note that a failure of a single pole to be attracted to a large cusp means that tip-splitting 
has occurred. This is the exact analog of tip-splitting in Laplacian growth. 

It is now time to relate the channel and radial geometries. We identify the typical scale 
in the radial geometry as C ^ W ^ Vq. On the one hand this leads to the scaling relation 
P — 1/x- On the other hand we use the result estabhshed in a channel, (4.26), with this 
identification of a scale, and find Tq — fo^- Comparing with (4.7) we find: 

This result leads us to expect two dynamical regimes for our problem. Starting from 
smooth initial conditions, in relatively short times the roughness exponent remains close to 
unity. This is mainly since the typical scale C is not relevant yet, and most of the poles that 
are generated by noise merge into a few larger cusps. In later times the roughening exponent 
settles at its asymptotic value, and all the asymptotic scaling relations used above become 
valid. We thus expect (3 to decrease from 1/(1 — /i) to an asymptotic value determined by 
X = l//? in (4.32): 

/3 = 1 + //Ri 1.35 ±0.03 . (4.33) 
The expected value of x is thus x — 0-74 ± 0.03. We tested these predictions in numerical 
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Figure 4.3: Panel a: a logarithmic plot of the correlation function F(y) of the interface versus 
the distance y between points. Panel b: Second derivative of the correlation function -F" (1) 
of the interface versus the distance 1 between points. 



simulations. We integrated Eq.(4.24), and in Fig. 4. 2 we display the results for the growth 
velocity as a function of time. After a limited domain of exponential growth we observe 
a continuous reduction of the time dependent exponent. In the initial region we get j3 = 
1.65 ±0.1 while in the final decade of the temporal range we find j3 = 1.35 ±0.1. We consider 
this a good agreement with (4.33). A second important test is provided by measuring the 
width of the system as a function of the radius, see Fig.4.2b. Again we observe a cross- 
over related to the initial dynamics; In the last temporal decade the exponent settles at 
X = 0.75 ± .1. We conclude that at times large enough to observe the asymptotics our 
predictions are verified. 

Finally, we stress some differences between radial and channel geometries. Fronts in 
a channel exhibit mainly one giant cusp which is only marginally disturbed by the small 
cusps that are introduced by noise. In the radial geometry, as can be concluded from the 
discussion above, there exist at any time cusps of all sizes from the smallest to the largest. 
This broad distribution of cusps (and scales) must infiucncc correlation function in ways that 
differ qualitatively from correlation functions computed in channel geometries. To make the 
point clear we exhibit in Fig.4.3a the structure function 

F{y) = ^{\R{x + y)-R{x)\^) (4.34) 

computed for a typical radial front, with x — R6. To stress the scaling region we exhibit the 
second derivative of this function in Fig.4.3b. The low end of the graph can be fitted well by a 
power law with a 0.6. This indicates that F{y) Ay + By"^^'^ . In a channel geometry 
we get entirely different structure functions that do not exhibit such scaling functions at all. 
The way to understand this behavior in the radial geometry is to consider a distribution 
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of cusps that remain distinct from each other but whose scales are distributed according to 
some distribution P'{t)- Let us define a distribution function P{t) ~ iP'{i) which give us 
the probabihty that a point on the circuit with the mean radius hes on the basis of the cusp 
with the size £. For each of these cusps there is a contribution to the correlation function of 
the form f{y,i) ~ ig{y/i) where g{x) is a scaling function, g{x) ~ a; for x < 1 and g{x) 
constant for x > 1. The total correlation can be estimated (when the poles are distinct) as 

F{y)^J2Pi^)^9{y/i) ■ (4.35) 
e 

The first derivative leaves us with J2iP{^)9'{y/^): and using the fact that g' vanishes for 
X > 1 we estimate F'{y) — Y^^yPi-^)- The second derivative yields F"{y) f» —P{y). Thus 
the structure function is determined by the scale distribution of cusps, and if the latter 
is a power law, this should be seen in the second derivative of F{y) as demonstrated in 
Fig. 4. 3. The conclusion of this analysis is that the radial case exhibits a scaling function that 
characterizes the distributions of cusps, P{£) ^ 

4.3 Conclusions 

The main purpose of this chapter was to find exponents of the problem and find connections 
between them. Using the main result of the channel case (dependence of the velocity on 
the ch annel size) we can find the acceleration of the fiame front in the radial case {Ro{t) = 
{const + tf) 

P^l + li . (4.36) 

where /i is the exponent for the dependence of the velocity on the channel size found 
before, /3 is the acceleration exponent and Rq is a mean radius of the fiame front. 

Dependence of the width of the flame front in the radial case W{t) on mean radius 
{Ro{t) = {const + tY) is 

W{t) ~ Ro{tY ■ (4.37) 

X = . (4.38) 

In summary, we demonstrated that it is possible to use information about noisy channel 
dynamics to predict nontrivial features of the radial evolution, such as the acceleration and 
roughening exponents. It would be worthwhile to examine similar ideas in the context of 
Laplacian growth patterns. 
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Chapter 5 
Laplacian Growth 



5.1 Introduction 

The problem of pattern formation is one of the most rapidly developing branches of nonlinear 
science today [1]. Of special interest is the study of the front dynamics between two phases 
(interface) that arises in a variety of nonequilibrium physical systems. If, as it usually hap- 
pens, the motion of the interface is slow in comparison with the processes that take place in 
the bulk of both phases (such as heat transfer, diffusion, etc.), the scalar field governing the 
evolution of the interface is a harmonic function. It is natural then, to call the whole process 
Laplacian growth. Depending on the system, this harmonic scalar field is a temperature 
(in the freezing of a liquid or Stefan problem), a concentration (in solidification from a su- 
persaturated solution), an electrostatic potential (in electrodeposition) , a pressure (in flows 
through porous media), a probability (in diffusion-limited aggregation), etc. 

The mathematical problem of Laplacian growth without surface tension exhibits a family 
of exact (analytic) solutions in terms of logarithmic poles in the complex plane. We show 
that this family of solutions has a remarkable property: generic initial conditions in channel 
geometry which begin with arbitrarily many features exhibit an inverse cascade into a single 
finger. 

In the absence of surface tension, whose effect is to stabilize the short-wavelength pertur- 
bations of the interface, the problem of 2D Laplacian growth is described as follows 

{dl + dl)u = . (5.1) 
u |r(t)= ,dnU |e= 1 ■ (5.2) 

Vn = dnU |r(t) . (5.3) 

Here u{x,y;t) is the scalar field mentioned, r{t) is the moving interface, E is a fixed 
external boundary, (9„ is a component of the gradient normal to the boundary (i.e. the 
normal derivative), and f„ is a normal component of the velocity of the front. 

We consider an infinitely long interface, obtained by a periodic continuation of the inter- 
face in the channel with periodic boundary conditions.Then We introduce a time-dependent 
conformal map / from the lower half of a "mathematical" plane, ^ = + irj, to the domain 
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of the physical plane, z = x + iy, where the Laplace equation 5.1 is defined as ^ — > z. We 
also require that f{t,^) ~ ^ for ^ — > ( — too. Thus the function z = f{t,() describes the 
moving interface. From Eqs. (5.1), (5.2), (5.3) for function f{t,^) we obtain the Laplacian 
Growth Equation 



^~di d^' ^ ' ^ ■ ^^■^> 

Now we will extend these results obtained for periodic boundary conditions to the more 
physical "no- flux" boundary conditions (no flux across the lateral boundaries of the channel) . 
This requires that the moving interface orthogonally intersects the walls of the channel. How- 
ever, unlike the case of periodic boundary conditions, the end points at the two boundaries 
do not necessarily have the same horizontal coordinate. This is also a periodic problem where 
the period equals twice the width of the channel. The analysis is the same as before, but 
now only half of the strip should be considered as the physical channel, whereas the second 
half is the unphysical mirror image. 

Let us look for a solution of Eq. (5.4) in the next form 

N 



fiU) = AC - ^r(^) - z^a,log(e'« - e^«'W), 

1=1 

N 

= 1 - A,-l < A < 1 , (5.5) 



1=1 



where r{t) is some real function of time, A is a real constant, ai is a complex constant, 
ii = Ci + iVi denotes the position of the pole with the number / and N is the number of poles. 
For the "no-flux" boundary condition we must add the condition that for every pole 
= + ^Vi with ai exists a pole — —Q + irji with aj. So for the function F{i^, t) — if{^, t) 



F{tU) = F{t^,t) (5.6) 
We want to prove that the flnal state will be only one flnger. 



5.2 Asymptotic behavior of the poles in the mathemat- 
ical plane 

The main purpose of this chapter is to investigate the asymptotic behavior of the poles in the 
mathematical plane. We want to demonstrate that for time t ^ oo, all poles go to a single 
point (or two points for no- flux boundary conditions). The equation for the interface is 

N 

/(e, i) = AC - iT{t) -iJ2ai log(e^« - e^^'W), 

1=1 

N 

^q;; = 1 - A,-l < A < 1 . (5.7) 

1=1 



65 



By substitution of Eq. (5.7) in the Laplacian Growth Equation 



dt 

we can find the equations of pole motion: 



^^y-^t — ) = 1 k=C-io , (5.8) 



and 



N 1 ^ 1 

const = T(t) + (1 - V a^) log — + V a;^log( a^) (5.9) 

fe=i t^i 



N N 

T = i- -^^a^Q;ilog(l-a^ai) + Co , (5.10) 

^ k=i 1=1 



where ai = e^^'. 

Prom eq. (5.9) we can find 

N N N 

Ci = (1 - A)t -^ailogai + ^ ^ o^o;; log(l -aj^ai) . (5.11) 

1=1 k=l 1=1 

Prom eqs. (5.10) and (5.11) we can obtain 

AT 

Im(^ ai log ai) — constant (5-12) 
1=1 

and 

1 + A 1 ^ 
t = (^^)r + -MH log «0 + Ci/2 , (5.13) 
^ ^ 1=1 

where ai is a constant, .^^(t) is the position of the poles, ai = e*^'^*^ , and is the portion of 
the channel occupied by the moving liquid. We will see that for r i— > oo we obtain one finger 
with wide ^ii. 

In Appendix A we will prove from eq.(5.10) that r i— > oo, if i i— > oo and if any finite time 
singularity does not exist. 

The equations of pole motion are next from eq. (5.9) 

canst = r + C + ^ a, log(l - e^^^'"^)), (5.14) 
const = Cfe + II(ariog I 1 - e'^^'^^^^ | +a'i arg(l - e^^^'^^^)), (5.15) 



I 



canst ^T + r)k + ^(aj log | 1 - e^(^'-^'=) | -a" arg(l - e^^^'"^"))), (5.16) 

Ci = Ci + iVhVi>0 . (5.17) 
ai — a'l + ia'l (5.18) 



66 



Let us examine 



arg(l - e^^^'-^*") = arg([l - e'^^'-^"-^ e'^^'^'^''^]) = arg[l - aike'"^'"] (5.19) 

Vik = Ci-Ck,aik = e-^'^^^+^'^^ (5.20) 
arg[l — a/fcC*''''*] is a singlevalued function of (pik, i.e. 

-| < arg[l - aike'^^^] < | . (5.21) 

From the eq. (5.16 ) the only way to compensate for the divergence of term r is that % i— > 
for T I— > oo, 1 < k < N. 

We want to investigate asymptotic behavior of poles r i— > cxo. To eliminate the divergent 
term log | 1 — e*^^*^"^*^) | we multiply eq. (5.16) by a'l and eq. (5.15) by a'j^ and take difference 

const = aiCk - a'Lr + - log I 1 - e'^^'"^^ | + 

l^k 

(«;«;, + aMarg(l-e'(«-«^))). (5.22) 

We have the divergent terms aj,'r in this equation. We may assume that for t i-^ oo, iV' 
groups of poles exist to eliminate the divergent terms {ipik ^ for all members of a group). 
Ni is the number of poles in each group, 1 < / < N' . For each group by summation of eqs. 
(5.22) over all group poles we obtain 

ccmst = af - af'V + Y.i{«' - log I 1 - e^^^'^^ | + 

l^k 

{ar'af + ar'af) arg(l - e^^^r-W))). (5.23) 

k 

c^r = T.^'k- (5.25) 

k 

We have no merging of these groups for large r and we investigate the motion of poles 
with this assumption 

I cr - cr i» vr + i<i,k<N. (5.26) 

For / k, r]f ^ 0, ipf^ = Cf - Ck we obtain 

-I gr 

log I 1 - e^^^r-^D \^ log I 1 - e^^^r-Cf) 1= log 2 + - logsin^ ^ (5.27) 

and 

gr 

arg(l - e^(r-C)) arg(l - 6^(^-0) = ^ + _ ^ . (5.28) 
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We choose n in Eq.(5.28) so that Eq.(5.21) is correct. 
Substituting these results to the eqs. (5.23) we obtain 



gr 

cu = «^"cr - <T + E[(«r«r - «) log i sin ^ i 

i^k ^ 

gr 

^{af ^ af" ar)^-\ (5.29) 

5.3 Theorem about coalescence of the poles 

From eqs. (5.29) we can conchidc 

(i) By summation of eqs. (5.29) (or exactly from eq. (5.12)) we obtain 

Y^a^'Ci'^ coast. (5.30) 

k 

gr 

(11) For I (^f^ |i— >• 0, 27r, we obtain log | sin ^ |i— > oo, meaning that the poles can not pass 
off each other; 

(iii) From (ii) we conclude that <| (/^f^ |< 27r 

(iv) From (i) and (iii), C,^^ i-^ oo is impossible; 

(v) In eq.(5.29) we must compensate the second divergent term. From (iv) and (iii) we 
can do it only if af" = oi^ = for all I. 

So from eq. (5.29) we obtain 







k 


(5.31) 




(5.32) 




(5.33) 




(5.34) 



For the asymptotic motion of poles in the group A^^ we obtain from eqs. (5.31), (5.32), 
(5.33), (5.34) taking leadind terms in eqs. (5.13), (5.14) 

2 

t , (5.35) 



A + 1 



= T + }_^a/-— -— — — . (5.36) 



y Vk + vi + KCk - 0) 

The solution to these equations is 



Vk^Vy^'^'\ (5.37) 
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^ik^^ik'e~-h'^''^-\ (5.38) 

4 = 0. (5.39) 
So we may conclude that for eliminating the divergent term we need 

Ni 

«r = E«^' = 0, (5.40) 

k 



af{l + X)>^ (5.41) 



for all I. 



5.4 Conclusions 

With the periodic boundary condition, eq.(5.40) is correct for all poles, so we obtain N' = 1, 
m = 1 and A^^ = N. 

Therefore the unique solution is 

Vk = vy^'' , (5.42) 

2 



Vik = Vik^e'T^)' , (5.43) 
4 = 0. (5.44) 

1 - > (5.45) 

With the no-flux boundary condition we have a pair of the poles whose condition of eq. 
(5.40) is correct so all these pairs must merge. Because of the symmetry of the problem these 
poles can merge only on the boundaries of the channel C = 0, ±7r. Therefore we obtain two 
groups of the poles on boundaries N' — 2, m — 1,2, Ni -\- — N , af ' + af^' = 1 — A. 

Consequently we obtain the solution (on two boundaries): 

^^-tf^'V^^-^*, (5.46) 

^S^ = ^..^^^'°e~+^^*, (5.47) 
C^^ = ; (5.48) 

1 

(2) (2),0 1+^ /r /|n^ 

Vk =Vk ^ ^ , (5.49) 
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Vlk = ^Ik^ ^' e "2 , (5.50) 

Cf ^ = ±^ ; (5.51) 

<'(! + A) > 0, (5.52) 

aT{l + A) > 0. (5.53) 



5.5 Appendix A 



We need to prove that r i— > cxd, if t i— > oo and if any finite time singularity does not exist. 
Tfiis is apparent if the second term in the next formula for r is greater than zero: 

N N 

T = t + [-- ^ ^aj^ai log(l - a^ai)] + Co , (5.54) 
^ fe=i 1=1 

where | a; |< 1 for all I. Let us prove it. 

N N ^ N N oo /— ^ \n 

-T^Y.Y.'^^i iog(i-zi^a;) = --^m^fe^'IK — hr~) 

^ k=l 1=1 ^ k=\ l=\ n=l 

-1 oo -1 AT N 

oE-(E^(^r)(E«K«/r) 

^ n=l fc=l Z=l 



-1 oo 1 JV TV 

^E-(E«K«0^)(E«K«0")>0 (5.55) 

^ n=l ^ ;=1 1=1 
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Chapter 6 
Summary 



The problem of flame propagation is studied as an example of unstable fronts that wrinkle 
on many scales. The analytic tool of pole expansion in the complex plane is employed to 
address the interaction of the unstable growth process with random initial conditions and 
perturbations. We argue that the effect of random noise is immense and that it can never 
be neglected in sufficiently large systems. We present simulations that lead to scaling laws 
for the velocity and acceleration of the front as a function of the system size and the level of 
noise, and analytic arguments that explain these results in terms of the noisy pole dynamics. 

We consider flame front propagation in channel geometries. The steady state solution 
in this problem is space dependent, and therefore the linear stability analysis is described 
by a partial integro-differential equation with a space dependent coefficient. Accordingly it 
involves complicated eigenfunctions. We show that the analysis can be performed to required 
detail using a flnite order dynamical system in terms of the dynamics of singularities in 
the complex plane, yielding detailed understanding of the physics of the eigenfunctions and 
eigenvalues. 

The roughening of expanding flame fronts by the accretion of cusp-like singularities is a 
fascinating example of the interplay between instability, noise and nonlinear dynamics that 
is reminiscent of self-fractalization in Laplacian growth patterns. The nonlinear integro- 
differential equation that describes the dynamics of expanding flame fronts is amenable to 
analytic investigations using pole decomposition. This powerful technique allows the devel- 
opment of a satisfactory understanding of the qualitative and some quantitative aspects of 
the complex geometry that develops in expanding flame fronts. 

Flame Propagation is used as a prototypical example of expanding fronts that wrinkle 
without limit in radial geometries but reach a simple shape in channel geometry. We show 
that the relevant scaling laws that govern the radial growth can be inferred once the sim- 
pler channel geometry is understood in detail. In radial geometries (in contrast to channel 
geometries) the effect of external noise is crucial in accelerating and wrinkling the fronts. 
Nevertheless, once the interrelations between system size, velocity of propagation and noise 
level are understood in channel geometry, the scaling laws for radial growth follow. 

The mathematical problem of Laplacian growth without surface tension exhibits a family 
of exact (analytic) solutions in terms of logarithmic poles in the complex plane. We show 
that this family of solutions has a remarkable property: generic initial conditions in channel 
geometry which begin with arbitrarily many features exhibit an inverse cascade into a single 
flnger. 
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